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A review of supcrsymmetry theory and phenomenology is presented. Topics dis- 
cussed include: gravity-mediated (SUGRA) and gauge-mediated (GMSB) super- 
symmetry breaking models; an overview of non-universal soft-supersymmetry- 
breaking masses and the resulting experimental implications; the phenomenology 
of and constraints on the possibility that a massive gluino is the lightest supersym- 
metric particle; current status of the phenomenology of the Higgs bosons of the 
minimal and next-to-minimal supersymmetric models; and, the signals for GMSB 
and/or R-parity- violating models coming from delayed decays of the lightest su- 
persymmetric partner of the Standard Model particles. 



1 Introduction 

There are good reasons why supersymmetry (SUSY) is regarded by many as 
the most attractive candidate for the next level of physics beyond the Standard 
Model (SM). First, there are important aesthetic motivations for SUSY, i) It 
is the only non-trivial extension of the Lorentz group at the heart of quantum 
field theory, ii) If SUSY is formulated as a local symmetry, it becomes a super- 
gravity (SUGRA) theory that reduces to general relativity in the appropriate 
limit, iii) SUSY appears in Superstrings. Second, although SUSY must be 
a broken symmetry since we have not yet observed supersymmetric particles, 
there are numerous indications that it will be realized at energies < few TeV. 
i) String theory solutions with a non-supersymmetric ground state are quite 
problematic, ii) SUSY solves the hierarchy problem, e.g. < 0- TeV) 2 

(as required to avoid a non-perturbative WW sector) via spin-1/2 loop can- 
cellation of spin-0 loop quadratic divergences provided the masses (generically 
denoted by tosusy) of supersymmetric partners are < few TeV. iii) The min- 
imal supersymmetric model (MSSM) leads to gauge coupling unification (at 
Mu ~ few x 10 16 GeV) if m SU SY £1-10 TeV. iv) In the MSSM, electroweak 
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symmetry breaking (EWSB) can be a natural result of renormalization group 
evolution (RGE) if m t ~ 175 GeV. v) In many cases, the lightest supersym- 
metric particle is a natural cold dark matter candidate. 

However, there are many phenomenologically viable supersymmetry mod- 
els. These and their experimental implications will be the primary subject of 
these lectures. The present review is an expanded version of Ref. 0. Full 
referencing for the more introductory material can be found there. Here, I will 
include detailed referencing only for some of the relatively recent developments 
and for the collider phenomenology discussions. 



2 The Minimal Supersymmetric Model 

The minimal supersymmetric extension of the Standard Model consists of tak- 
ing the Standard Model and adding the corresponding supersymmetric part- 
ners. The SM particles and corresponding sparticle mass eigenstates are listed 
below. 

[u,d,c,s,t,b] g W ± ,H ± y,Z,H°,H° 

u, d, c, s, t, b 
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In addition, the MSSM contains two hypercharge Y = ±1 Higgs doublets, 
which is the minimal structure for the Higgs sector of an anomaly-free super- 
symmetric extension of the Standard Model. The supersymmetric structure 
of the theory also requires (at least) two Higgs doublets to generate mass for 
both "up" -type and "down" -type quarks (and charged leptons). The physical 
Higgs mass eigenstates then comprise two CP-even Higgs bosons, h°, H° (with 
m/jO < m H o), a CP-odd Higgs boson, ^4°, and a charged Higgs pair, i? . 

In the MSSM, all renormalizable supersymmetric interactions consistent 
with (global) B—L conservation are included. They are parameterized by: 

• SU(3)xSU(2)xU(l) gauge couplings, g 3 , g 2 = g, and g x = \J\g'\ 

• (complex) Higgs-quark Yukawa coupling matrices, y u , yd, and y e ; and 

• a supersymmetric Higgs mass parameter, [i. 

Since the MSSM is a model of three generations of quarks, leptons and their su- 
perpartners, y u , yd, and y e are complex 3x3 matrices. However, as discussed 
below, not all these degrees of freedom are physical. 

As a consequence of B—L invariance, the MSSM possesses a multiplicative 
R-parity invariance, such that all the ordinary Standard Model particles have 



even R-parity, whereas the corresponding supersymmetric partners have odd 
R-parity. 

In the MSSM, supersymmetry breaking is accomplished by including the 
most general renormalizable soft-supersymmetry-breaking terms consistent with 
the SU(3)xSU(2)xU(l) gauge symmetry and R-parity conservation. These 
terms parameterize our ignorance of the fundamental mechanism of super- 
symmetry breaking. It is here where most of the new supersymmetric model 
parameters reside. As already noted, if supersymmetry is relevant for explain- 
ing the scale of electroweak interactions, then the mass parameters introduced 
by the soft-supersymmetry-breaking terms should be in the TeV range or be- 
low. 

The soft-supersymmetry-breaking terms of the MSSM include: 

• gaugino Majorana masses M 3 , M 2 , and Mi, associated with the SU(3), 
SU(2), and U(l) subgroups of the Standard Model; 

• diagonal Higgs squared-mass terms, tti 2 h u1 m ff d ' 

• an off-diagonal Higgs squared-mass term, b (the "B"-term); 

• (hermitian) squark and slepton squared-mass matrices, niq, m^, m|, 
ml, m|; 

• (complex) matrix ^4-parameters, a u , ad, a e , parameterizing Higgs-squark- 
squark and Higgs-slepton-slepton trilinear interactions. 

As above, the quantities specified by boldface letters are 3x3 matrices (in 
generation space). Note that after minimizing the Higgs potential, one can 
trade in the three Higgs parameters mjj , and b for two neutral Higgs 

field vacuum expectation values v u and Vd and one of the physical Higgs masses. 
Equivalently, one fixes v 2 = v 2 . + v 2 d = Am^/g 2 = (246 GeV) 2 and introduces 
the parameter tan/3 = v u /vd- 

There are many unphysical degrees of freedom among the MSSM pa- 
rameters listed above. These can be eliminated by expressing interaction 
eigenstates in terms of the mass eigenstates, with an appropriate redefini- 
tion of the MSSM fields to remove unphysical phases. In the end, the MSSM 
possesses 124 truly independent parameters. Of these, 18 parameters cor- 
respond to Standard Model parameters (including the QCD vacuum angle 
$QCd), one corresponds to one of the Higgs sector masses (the analogue of 
the Standard Model Higgs mass), and 105 are genuinely new parameters of 
the model. These include: five real parameters and three CP-violating phases 
in the gauge/gaugino/Higgs/higgsino sector, 21 squark and slepton masses, 



36 new real mixing angles to define the squark and slepton mass eigenstates 
and 40 new CP-violating phases that can appear in squark and slepton in- 
teractions. The most general R-parity conserving minimal supersymmetric 
extension of the Standard Model (without additional theoretical assumptions) 
will be denoted henceforth as MSSM-124. 

MSSM-124 is phenomenologically viable only for very special domains in 
the full parameter space such that the following possible problems are absent: 
(i) non-conservation of the separate lepton numbers L e , L M , and L T ; (ii) un- 
supprcsscd flavor-changing neutral currents (FCNC's); and (iii) new sources of 
CP-violation that are inconsistent with the experimental bounds. 

There are a number of approaches for reducing the parameter freedom 
of MSSM-124 in such a way as to avoid these problems. There are two low- 
energy approaches. First, one can assume that the squark and slepton squared- 
mass matrices and the matrix ^4-parameters are proportional to the 3x3 unit 
matrix (horizontal universality). Alternatively, one can simply require that 
all the aforementioned matrices are flavor-diagonal in a basis where the quark 
and lepton mass matrices are diagonal (flavor alignment). In cither case, L e , 
and L T are separately conserved, while tree- level FCNC's are automatically 
absent. In both cases, the number of free parameters characterizing the MSSM 
is substantially less than 124. Both scenarios are phenomenologically viable, 
but there is no strong theoretical basis for either scenario. 

In high-energy approaches, one treats the parameters of the MSSM as run- 
ning parameters and imposes a particular structure on the soft-supersymmetry- 
breaking terms at a common high-energy scale [such as the GUT scale Mjj or 
the reduced Planck scale M P = (87rGW)~ 1/2 _ 2 .4 x 10 18 GeV]. Using the 
renormalization group equations, one can then derive the low-energy MSSM 
parameters. The initial conditions (at the appropriate high-energy scale) for 
the renormalization group equations depend on the mechanism by which super- 
symmetry breaking is communicated to the effective low-energy theory. One 
bonus of such an approach is that one of the diagonal Higgs squared-mass pa- 
rameters is typically driven negative by renormalization group evolution. Thus, 
electroweak symmetry breaking is generated radiatively, and the resulting elec- 
troweak symmetry-breaking scale is intimately tied to the scale of low-energy 
supersymmetry breaking. 

In the context of the MSSM, the high-energy approach is especially well- 
motivated given the fact that the MSSM particle content yields relatively exact 
gauge coupling unification (for sparticle masses < 1 — 10 TeV). If the gauge 
couplings become universal at Mjj, then it is quite natural to suppose that 
the soft supersymmetry breaking parameters could also take a simple form 
at Mjj. A number of points concerning gauge coupling unification should be 



emphasized, a) Pcrturbative unification requires that the number of families 
in the < 1 TeV mass range be < 4. b) Unification occurs only if there are 
exactly 2 Higgs doublets, as in the MSSM, and possibly additional singlets, 
below Mij. Coupling unification in the presence of additional doublets, triplets 
and/or higher Higgs representations is only possible if appropriate additional 
matter fields with mass above 1 — 10 TeV and below Mjj are introduced, c) The 
currently preferred smaller a s (mz) values are most consistent with coupling 
unification if the 'average' sparticle mass msuSY is above 1 TeV. d) The fact 
that the predicted MSSM value of Mjj is substantially below Mp may prove 
difficult to accommodate in specific string models. 

To summarize, the word "minimal" in MSSM refers to the minimal particle 
spectrum and the associated R-parity invariance. The MSSM particle content 
must be supplemented by assumptions about the origin of supersymmetry- 
breaking that lie outside the low-energy domain of the model. The corre- 
sponding parameter freedom is bad if you want to know ahead of time exactly 
how to discover supersymmetry. Machine builders and experimentalists must 
be prepared for a wide range of possibilities as to how to see and fully explore 
SUSY. But, it is good in the sense that once the 105 new parameters of the 
MSSM have been experimentally determined you will have learned a great deal 
about the model. For example, in models where supersymmetry-breaking is 
mediated by gravitational interactions one can hope to discover an underlying 
organization for the 105 parameters (by evolving the \ow-E parameters up to 
Mjj or Mp) that will define the GUT/String model visible sector properties. 
Different models lead to vastly different \ow-E parameters and expectations. 
The rest of the talk reviews some of the many possibilities discussed to date. 

3 Constrained MSSM's: mSUGRA and GMSB 

In these theories, one posits a theory of dynamical supersymmetry-breaking 
(DSB), where the DSB occurs in a sector that is distinct from the fields of 
the MSSM. The supersymmetry-breaking inherent in the DSB sector is subse- 
quently transmitted to the MSSM spectrum by some mechanism. 

Two theoretical scenarios have been examined in detail: gravity-mediated 
and gauge-mediated supersymmetry breaking. In these approaches, supersym- 
metry is spontaneously broken, in which case a massless Goldstone fermion, 
the goldstino, arises. Its coupling to a particle and its superpartner is fixed by 
the supersymmetric Goldberger-Trieman relation: 

C ia t = -j.f a d^G a + h.c, (1) 
where is the supercurrcnt, which depends bilinearly on all the fermion 



boson superpartner pairs of the theory and G a is the spin- 1/2 goldstino field 
(with spinor index a). yF is the scale at which supersymmetry-breaking 
occurs in the hidden sector (typically, y/F ^S> mz)- When gravitational effects 
are included, the goldstino is "absorbed" by the gravitino (93/2), the spin- 
3/2 partner of the graviton. By this super-Higgs mechanism, the goldstino is 
removed from the physical spectrum and the gravitino acquires a mass (m 3 / 2 ). 
In models where the gravitino mass is generated at tree-level one finds 

m3/2 = v£p- (2) 

The helicity ±i components of the gravitino behave approximately like the 
goldstino, whose couplings to particles and their superpartners are determined 
by Eq. (jl]). In contrast, the helicity ±| components of the gravitino couple 
with gravitational strength to particles and their superpartners, and thus can 
be neglected. 

In many models, the DSB sector is comprised of fields that are completely 
neutral with respect to the Standard Model gauge group. In such cases, the 
DSB sector is also called the "hidden sector" . The fields of the MSSM are 
said to reside in the "visible sector", and the model is constructed such that 
no renormalizable tree-level interactions exist between fields of the visible and 
hidden sectors. A third sector, the so-called "messenger sector", is often em- 
ployed in models to transmit the supersymmetry-breaking from the hidden 
sector to the visible sector. However, it is also possible to construct models in 
which the DSB sector is not strictly hidden and contains fields that are charged 
with respect to the Standard Model gauge group. 

3.1 Gravity-mediated super symmetry breaking 

All particles feel the gravitational force. In particular, particles of the hidden 
sector and the visible sector can interact via the exchange of gravitons. Thus, 
supergravity (SUGRA) models provide a natural mechanism for transmitting 
the supersymmetry-breaking of the hidden sector to the particle spectrum of 
the MSSM. In models of gravity-mediated supersymmetry breaking, gravity is 
the messenger of supersymmetry-breaking. The resulting low-energy effective 
theory below the Planck scale consists of the unbroken MSSM plus all possible 
soft-supersymmetry-breaking terms. 

For example, the gaugino masses could arise primarily from an effective 
Lagrangian term of the form 

£= [ d 2 9W a W b ^ + h.c. ~ { -^\ a \ b , (3) 



where the A 0; h (a, b — 1,2, 3) are the gaugino fields and F$ is the auxiliary com- 
ponent of a hidden sector superfield $ that appears linearly in the gauge kinetic 
function. The scale Mp appears in the denominator when communication be- 
tween the hidden sector and the MSSM sector is via gravitational interactions. 
F in the previous section should be identified with (F$). From Eq. (||) we see 
that in order for gauginos to have mass 0(mw), VF ~ V m wMp ~ 10 10 GeV 
is required. Such a large \f~F value implies that the gravitino mass, Eq. (|J), is 
larger than the mass of the lightest neutralino, Xi- Since the gravitino is also 
very weakly interacting when \J~F~ is large, Eq. ([[J), it becomes phenomenolog- 
ically irrelevant. 

In a minimal supergravity (mSUGRA) framework, the soft-supersymmetry 
breaking parameters at the Planck or GUT scale take a particularly simple 
form. First, the gaugino mass parameters as well as the gauge couplings are 
presumed to unify at some scale Mx (= Mjj or Mp?): 

M x (M x ) = M 2 (M x ) = M 3 (M x ) = m 1/2 . (4) 

[For example, in the context of Eq. (0), the field $ of the hidden sector would 
be presumed to be a flavor singlet, implying oc c5 a b] Eq. (|J) implies that 
the low-energy gaugino mass parameters satisfy: 

M 3 = %M % ~ 3.5M 2 , Mi = | tan 2 6 W M 2 ~ 0.5M 2 . (5) 
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Second, the scalar squared masses and the yl-parameters are taken to be flavor- 
diagonal and universal: 

(M P ) = m=L(M P ) = m|(M P ) = m£ (M P ) = m|(M P ) = m 2 l , 
™ 2 hAM p ) = m 2 Hd (M P ) = m 2 , 

a f (M P ) = A y f (M P ) , f = u,d,e, (6) 

where 1 is the 3x3 identity matrix in generation space, and the yf are the 
Higgs-fermion Yukawa coupling matrices. Finally, /i, b and the gaugino mass 
parameters are assumed (rather arbitrarily) to be initially real at the high 
scale. 

It is easy to count the number of free parameters of mSUGRA. The low- 
energy values of the MSSM parameters are determined by the MSSM renor- 
malization group equations and the above initial conditions [Eqs. (^) and (H)]. 
Thus, the mSUGRA model is fixed by the following parameters: 18 Stan- 
dard Model parameters (excluding the Higgs mass), mo, to^, Aq, tan/3, and 
sgn(/i) for a total of 23 parameters. This result was obtained by trading in 



the parameters b and /j, 2 for v 2 (which is counted as one of the 18 Standard 
Model parameters) and tan/3. In this procedure the sign of /i is not fixed, 
and so it remains an independent degree of freedom. The requirement of ra- 
diative electroweak symmetry breaking imposes an additional constraint on 
the possible range of the mSUGRA parameters. In particular, one finds that 
1 < tan/3 < rut/nib- In principle, the mass of the gravitino, m 3 / 2 (or equiva- 
lently, \AF), is an additional independent parameter. But, as discussed above, 
the gravitino does not play any role in mSUGRA phenomenology. 

Clearly, mSUGRA is much more predictive than MSSM-124. In particu- 
lar, one has only four genuinely new parameters beyond the Standard Model 
(plus a two- fold ambiguity corresponding to the sign of /i), from which one can 
predict the entire MSSM spectrum and its interactions. It will be convenient 
to reference a number of specific mSUGRA parameter sets. The first is the 
'Standard' Snowmass96 Point, defined by too = 200 GeV, m 1 / 2 — 100 GeV, 
Aq = 0, tan/? = 2, fi < 0. Other interesting parameter sets include the 
String/SUGRA motivated No-Scale (toi/2 7^ 0, too = Aq = 0) and Dilaton 
or Dilaton-Like (to!/ 2 = — Aq = \/3mo) models. The latter are perhaps par- 
ticularly worthy of considering since to deviate significantly from dilaton-likc 
boundary conditions in Calabi-Yau and Orbifold models requires that super- 
symmetry breaking be overwhelmingly dominated by moduli fields other than 
the dilaton (sin go idstmo < 1). This is because in the absence of loop correc- 
tions the dilaton boundary conditions apply with m 1 / 2 — to 3 / 2 sin go idstino and 
it is only when m 3 / 2 sin # go idstino is small compared to loop corrections that we 
can move away from dilaton boundary conditions. Vacuum stability against 
color breaking, might be a problem for this model, but it currently appears 
that metastability and/or non-perturbative effects can cure the apparent diffi- 
culty. In both the no-scale and dilaton cases, two parameters (usually chosen 
to be TO!/ 2 and tan/3) and the sign of fj, are sufficient to completely specify the 
model. 

The mSUGRA initial conditions of Eq. (||) correspond to a minimal SUGRA 
framework (specifically, the kinetic energy terms for the gauge and matter 
fields are assumed to take a minimal canonical form). Since there is no theo- 
retical principle that enforces such a minimal structure, perturbations of the 
mSUGRA initial conditions consistent with phenomenology (e. g, avoiding the 
generation of dangerous FCNCs) have been considered For example, 

one can introduce separate mass scales for the Higgs and squark/slepton soft- 
supersymmetry-breaking masses. One can also introduce non-universal scalar 
masses, but restrict the size of the non-universal terms to be consistent with 
phenomenology Finally, gaugino mass unification is not required by either 
theory or phenomenological constraints. Specific examples of all these devia- 



tions arise in string- inspired models Q. Models which violate the mSUGRA 
initial conditions that appear to be of particular interest include: 

• Strange-event-motivated models. For example, the phenomenology of 
Mi ~ M 2 [in contrast to M x ~ 0.5M 2 predicted by Eq. f])] has recently 
been advocated || in order to provide a possible explanation for the 
famous CDF ee77 event. 

• Scenarios with M% < Mi , resulting in the Xi an d Xi both being wino-like 
and nearly degenerate in mass. 

• The light gluino scenario, in which the Mi are set to zero. 

• Scenarios in which non-universality for the M, is such that the gluino 
has substantial mass, but is the lightest of the gauginos. 

Non-Universal Gaugino Masses 

Non-universal gaugino masses are particularly interesting from a phenomeno- 
logical viewpoint. Models with substantial theoretical motivation include the 
following. 

• The O-II Orbifold String model with all matter fields in the n = — 1 
untwisted sector, and supersymmetry breaking dominated by the overall 
size modulus (as opposed to the dilaton). The phenomenology of this 
model is explored in Ref. 

This is the only known string-based model in which the limit of pure modulus 
(no dilaton) supersymmetry breaking is possible without charge and/or color 
breaking. The gaugino masses arise from 1-loop terms generic to any theory; 
the resulting Mjj boundary conditions take the form 

M%~V3m 3/2 [-(b a + 5 GS )Kr)], m 2 Q = mj /2 [-6 GS K'}, A = 0, (7) 

where the b a are the gauge RGE coefficients (63 = 3, b 2 = — 1, &i = —33/5), Sqs 
is the Green-Schwarz mixing parameter (a negative integer in the O-II model, 
with (5gs = ~3, — 4, — 5 preferred), and 77 = ±1. The one-loop estimates of 
K = 4.6 x 10 -4 and K' = 10~ 3 imply mj, m~3> M a . We emphasize that the 
above relation of M a to b a , Sqs that arises at one-loop is more general than the 
O-II model; in particular, it is likely to survive non-perturbative corrections. 
However, the relation between K and K' is model-dependent. 

• F-term breaking with F 7^ SU(5) singlet. The phenomenology of these 
models is discussed in Ref. m. 
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Table 1: M a at Afy and mz for the four F$ irreducible representations and in the O-II 
model with Sqs ~ —4. 



In these models, it is supposed that the hidden sector is not completely hidden; 
it contains some fields with ordinary quantum numbers. In particular, referring 
to Eq. (Q), suppose that 4> is not a singlet field. In general, it is only required 
that 

$, F$ £ (24x24) symmctri c = 1 © 24 © 75 200 . (8) 

For higher representations, only F$ components that are 'neutral' with respect 
to the SU(3), SU(2), U(l) gauge groups can acquire a vev if these groups are to 
remain unbroken after supersymmetry breaking. The result is that (i 7 $) a ( ) = 
CaSab, with c a depending on the representation (an arbitrary superposition of 
representations also being possible). The relative gaugino masses are tabulated 
in Tabic g. 

To obtain a first appreciation for the phenomenological implications we 
give in Table the low-energy gaugino masses for these cases, assuming the 
Snowmass96 point m~ ~ 285 GeV value for the gluino mass. The fact that 
to~ ~ min(Mi,M2) and m~± ~ M 2 (after mass matrix diagonalization) im- 
plies that in the 200 and O-II models (which have M 2 < Mi) m~± ~ m~ 

(and both are winos). Of course, the relative masses depend upon tan/3 as 
illustrated in Ref. S. In the O-II model, the dependence of the ino masses 
on Sqs, illustrated in Fig. |l|, is also quite illuminating. Note that for a nar- 
row range of 5gs (near the theoretically preferred values), the gluino will be 
the LSP! The detailed phenomenological implications of these results will be 
outlined later. 

• The light-gluino scenario. 

In this scenario 0, the soft-supersymmetry-breaking gaugino masses Mi are 
taken to be zero at tree-level and only arise at one-loop. The result is that Xi ~ 





1 


24 


75 


200 


O-II 
Sgs = —4.7 


772 — 
9 

ra~ 

ra- 
ti 

ra- 
ti 

m~ 

b L 

ra~ 

ra- 
in. 

ra- 
th 

ra~ 


285 
302 
255 
315 
266 
303 
207 
216 
44.5 


285 
301 
257 
321 
276 
303 
204 
229 
12.2 


287 
326 
235 
351 
307 
309 
280 
305 
189 


288 
394 
292 
325 
264 
328 
437 
313 
174.17 


313 
big 
big 
big 
big 
big 
big 
big 
303.09 


ra~ 


97.0 


93.6 


235 


298 


337 


ra— + 

Xl 


96.4 


90.0 


240 


174.57 


303.33 


m~+ 
x 2 


275 


283 


291 


311 


big 


ra h o 


67 


67 


68 


70 


82 



Table 2: Sparticle masses for the Snowmass96 comparison point. 



7 with m~ ~ 0.1 — 1.5 GcV (from one-loop squark-quark and wmo/higgsino- 
Higgs/ vector boson diagrams) and m~ ~ 10 — 600 MeV (from one-loop squark- 
quark diagrams). The lightest supersymmetric particle will be the lightest 
'R-hadron', which is likely to be the 7?° = gg with mass of order 1.5 GeV. 
It is claimed || that this model is highly disfavored on the basis of 2j and 
4j event shapes in Z decay, which, for example, constrain the running of a s 
via b = (11/6)C A /C F - {2/3)n f T F /C F . A light g acts like 3 new flavors, 
increasing rif from ~ 5 to ~ 8. Another immutable prediction of the model is 
that m~± < mw- Thus, LEP-2 running at yfs ~ 190 GeV should soon find 

Xi 

evidence for the chargino (despite the possibly complicated decay patterns) or 
the model will be ruled out. It is important to note that these searches do not 
constrain the possibility noted in the previous item of a heavy gluino-LSP. 

3.2 Gauge-mediated supersymmetry-breaking 

The theory of gauge-mediated supersymmetry breaking (GMSB) || posits 
that supersymmetry breaking is transmitted to the MSSM via gauge forces. 
Two model-building approaches have been recently explored in the literature. 
First, in hidden-sector models, the GMSB model consists of three distinct sec- 
tors: a hidden (DSB) sector where supersymmetry is broken, a "messenger" 




sector consisting of messenger fields with SU(3) x SU(2) x U(l ) quantum num- 
bers, and a sector containing the fields of the MSSM |L0|. The direct coupling 
of the messengers to the hidden sector generates a supersymmetry-breaking 
spectrum in the messenger sector. Second, in models of "direct gauge medi- 
ation" (sometimes called direct-transmission models) jy], the GMSB model 
consists only of two distinct sectors: the DSB sector (which also contains the 
messenger fields) and the sector of MSSM fields. In both classes of models, 
supersymmetry-breaking is transmitted to the MSSM via SU(3)xSU(2)xU(l) 
gauge interactions between messenger fields and the MSSM fields. In par- 
ticular, supersymmetry-breaking masses for the gauginos and squared-masses 
for the squarks and sleptons arise, respectively, from one-loop and two-loop 
diagrams involving the virtual exchange of messenger fields. 

In simple hidden-sector GMSB models, gaugino masses have the same 
relative values as if they were unified as specified in Eq. ([|) (with Mx of 
order My) even though the actual initial conditions are set at scale M m . For 
example, in the minimal model where supersymmetry-breaking effects in the 
messenger sector are small compared to M m , one finds the following result for 



the gaugino Majorana mass terms [ |To| : 



Mi(M m ) = fc^„ lff ^^A, (9) 

where &2 = &3 = 1, k\ = 5/3, and N m denotes the effective number of messen- 
gers in the messenger sector — N m < 4 is required to avoid Landau poles. A 
is an effective mass scale that depends on the details of the GMSB model. In 
addition, scalar masses are expected to be flavor-independent since the medi- 
ating gauge forces are flavor-blind. In the simplest hidden-sector models they 
take the form: 

ml{M m ) = 2A»JW (^) |> (^) 2 (10) 

with c 3 = 4/3 (triplets) C2 = 3/4 (weak doublets), and Ci = § (y) 2 in the 
normalization where Y/2 = Q — T3. Thus, degeneracy among families is broken 
only by effects of order quark or lepton Yukawa couplings. In Eqs. (||) and ([l0|), 
M m / A > 1 is required to avoid negative mass-squared for bosonic members of 
the messenger sector; M m /A > 1.1 is preferred to avoid fine-tuning, for which 
f(A/M m ) ~ 1 and 1 < g(A/M m ) < 1.23. For M m /A > 2, 1 < g(A/M m ) < 
1.045. In the approximation g = f = 1, one finds (at tree-level) m~ : m~ : 

mj : Mi — 11.6 : 2.5 : 1.1 : >/N m , ^ which implies that the next to lightest 
supersymmetric particle (the goldstino being the lightest) is the Xi — B l° r 
N m — 1 or the i^s for N m > 2. Note that supersymmetric scalar and gaugino 
states with the same color quantum numbers have masses of roughly the same 
magnitude. 

In order that all superpartner masses be < 1 TeV, it follows [e.g., from 
Eq. @] that A < 100 TeV. For example, for N m = 1 

where Mi is the low-energy value of the U(l)-gaugino soft-supersymmetry- 
breaking mass and the gluino mass is m~ ~ 7M%. The masses of various 
sparticles in the case of N m — 1 are plotted as a function of A (taking / = g = 
1) in Fig. |. 

^Evolution corrections ]l2| in going from scale M m to the <^ 1 TeV scale do not destroy this 
ordering unless M m is very large, although they do affect the precise numerical ratios. 
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Figure 2: m~ n , m~ + , m~ , m~ , m~, m~ and m~ as functions of A in the N m = 1 GMSB 
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scenario, taking / = g = 1 and M m = 1.1A. 



As for the other parameters of low-energy supersymmetry, the A-parameters 
are suppressed, while the generation of fi and b is quite model-dependent (and 
lies somewhat outside the standard ansatz of gauge-mediated supersymmetry 
breaking). The initial conditions for the soft-supersymmetry- breaking running 
parameters are fixed at M mi which characterizes the average mass of messen- 
ger particles. In principle, M m can lie anywhere between (roughly) A and 
10 16 GeV (in models with larger values of M m , supergravity- mediated effects 
would dominate the gauge-mediated effects). Thus, the minimal GMSB model 

contains 18 Standard Model parameters, A, tan (5 and arg(/i) [after trading 
in b and \n\ 2 for v and tan/3]. There is also a weak logarithmic dependence 
on M m , which enters through renormalization group running. We thus end up 
with 22 free parameters. Clearly, the minimal GMSB approach is even more 
predictive than the mSUGRA model. 

When the GMSB approach is extended to incorporate gravitational phe- 
nomena, supergravity effects will also contribute to supersymmetry breaking. 
However, in models of gauge-mediated supersymmetry breaking, one usually 
chooses the model parameters in such a way that the virtual messenger ex- 
change diagrams dominate the effects of the direct gravitational interactions 
between the hidden and visible sectors. However, of critical importance is the 
fact that the super-Higgs effect becomes operative, and the gravitino absorbs 
the goldstino and becomes massive. Using Eq. (|h, it is convenient to rewrite 



the formula for the gravitino mass as follows: 



The scale of supersymmetry breaking, v F, is a very crucial parameter in 
GMSB models. It is highly model-dependent. In hidden-sector models, val- 
ues of y/F > 10 3 -10 4 TeV are required in order to generate sufficiently large 
supersymmetry-breaking in the sector of MSSM fields JlJ,^5 16 . In particular, 
one can derive ]l4||ll| the approximate inequality 

T ^>,(^)('), (13) 



V 2000 TeV I ~ J \100 GeV J \27.5 

where / ~ (if^) ~ 2-5 x 10 4 /g^, with g m being the coupling of the gauge 

group responsible for communicating (via two-loop diagrams) supersymmetry 
breaking from the hidden sector to the messenger sector. Pcrturbativity would 
require g m < 1. For to- > 45 GeV (the rough LEP limit) and / = 2.5 x 10 4 , 

VF > 5200 TeV (VF > 2600 TeV) for N m = 1 (N m = 4). g In direct- 
transmission models, the two-loop communication between the hidden and 
messenger sectors is eliminated, / in Eq. (|k|) is effectively of order unity, and 
y/~F can be as low as 100 TeV in phenomenologically viable models. 

Eq. ( |l2| ) implies that the gravitino can be quite light for moderate y/~F 
values > 100 TeV. For VF values such that to— is smaller than to~ , m~ , the 

— 93/2 Xi' Ir' 

<?3/2 w ih be the LSP. An upper bound on \/~F then arises in R-parity conserving 
GMSB models from the fact that the gravitino LSP is stable, and thus is 
also a candidate for dark matter. While very light gravitinos (eV masses) 
will not contribute significantly to the total mass density of the universe, if 
the gravitino LSP is too heavy (m~ ^ > few keV) its relic mass density in 

most early universe scenarios will be greater than the critical density p7| , p^| 
in conflict with observation. ^ 

From Eq. ((l^), this means that y/~F values above ~ 3000 TeV are dis- 
favored. As previously noted, the helicity ±i components of g 3 / 2 behave 
approximately like the goldstino. In contrast to SUGRA models, F^ 1 is 

c One should allow for a factor of two or three uncertainty in the approximate lower bound. 
^Although the LSP is an ideal candidate for being a major component of the dark matter, one 
need not use this as a constraint to restrict y/F. It may turn out that the main component of 
the dark matter has another source. Some examples are: the QCD axion, its supersymmetric 
partner (the axino) or the lightest stable particle in the GMSB messenger sector. 



now significantly larger, so that the goldstino coupling to the particles of the 
MSSM [Eq. (|l])] is significantly enhanced relative to its SUGRA-model cou- 
pling strength. Thus, in GMSB models, the gravitino is very likely to be the 
LSP and will play a phenomenologically crucial role. 

However, minimal GMSB is not a fully realized model. The sector of 
supersymmetry-breaking dynamics can be very complex, and it is fair to say 
that no complete model of gauge-mediated supersymmctry yet exists that is 
both simple and compelling. In particular, the various direct transmission 
models are very different from one another and from the simple hidden-sector 
model. Relations between sparticle and gravitino masses vary tremendously. 
As an extreme example, in the model of Ref. |Q the gluino is the LSP and 
the gravitino is sufficiently heavy as to be phenomenologically irrelevant. 



4 Constraints on the gluino-LSP possibility 

We have seen that it is possible in both SUGRA and GMSB models that the 
lightest R-hadron containing the (massive) gluino is the LSP. Thus, this rather 
unusual possibility deserves serious consideration. 

First, there are many limits on heavy stable charged particles based on 
mass spectrometer searches for heavy isotopes of hydrogen and oxygen and 
accelerator fixed target searches for heavy new charged particles. Thus, we 
assume that it is the neutral Ro — gg that is the LSP. We assume that other 
neutral states such as the p° = g(uu + dd) and charged states such as the 
p + = gud (and baryonic equivalents) have short lifetimes for decay to the R°. 

Limits on a strongly interacting neutral LSP have been considered in 
Refs. |l9| , ^0| , |l8| . There are many differences in these three treatments. The 
assessment of constraints is a very messy subject and is in a state of flux. 
Nonetheless, the topic is sufficiently interesting to warrant a brief overview. 

As a first ingredient, we need an estimate for the relic density of R°'s. 
Ref. pp| , using standard freeze-out arguments and <r ann ~ a 2 s /m 2 R0 , finds 
^2- - 10~ 3 m R o( TeV), whereas Ref. pi obtains ^ _ 10~ 7 m R a( TeV). 
A factor of 10 3 in this difference appears to be due to the much larger an- 
nihilation cross section [<7 an n ~ 40 mb x (m pro ton/wflo) ] implicitly assumed 



in Ref. 18 . Even for the larger relic density of Ref. p3|, one finds 



PB 

~ ( Yt§v ) ' hnplying that the R° cannot be the source of dark matter 
if m^jo < 1 TeV. This result, in turn, implies that if m R o < 1 TeV then it is 
unlikely that the galactic halo is entirely made of i?°'s. 

Given a density of i?°'s, an important ingredient in obtaining limits is 
the flux in the vicinity of the earth of the i?°'s. If the halo density (p — 
0.4 GeV/cm 3 ) is dominated by i?°'s, implying n R o = p/m^o, and using 



the velocity of movement of the earth through the halo (v R o ~ 2 — 3 x 
10 7 cm/s), the flux of R 0, s from the galactic halo would be cf> R o = n R ov R a ~ 

10 4 ( 1 TcV ) cm' 2 s -1 . On the other hand, if we employ the relic density 
^ ~ 10~ 3 m fi o( TeV) from Rcf. §0| with n B ~ 2 x KT 7 cm" 3 , one finds 
(using the same v R o) <j) R o ~ 6 x ( t j^v ) cm^s -1 , *- e - much smaller for 
m R o ^ 1 TeV. 

Other important ingredients in obtaining limits on the R° are the R° cross 
sections. Ref. [§0) estimates a{R°p) ~ 10 mb and a ann (R°R ) ~ a 2 Jm R0 (a 
larger value applying if the kinetic energies of the i?°'s are below a GeV or so). 
Ref. 1 18 suggests that the leading contribution to a(R°p) would come from 
glueball exchange: <r(R°p) ~ 40(rn 7r /m g i uc ban) 2 mb. Since the lightest glueball 
is ~ 10 times heavier than a pion, a{R°p) would then be < 0.4 mb. 

Limits on the R° derive from: (1) dark-matter (wimp) balloon and un- 
derground mine detection experiments; (2) capture and co-annihilation in the 
earth Q (to produce j/'s detected in proton decay detectors, at least as sec- 
ondary products of primary annihilation) or capture and co-annihilation in 
the halo |^0| to produce energetic 7's that are detected; (3) limits on heavy 
isotopes created by the flux of i?°'s impacting the earth. We consider each in 
turn. 

(1) Ref. (2^] considers only underground wimp detectors. First, they argue 
that since the i?°-nuclear cross section is > 10 10 times larger than the generic 
heavy neutrino cross section, the wimp detectors would certainly have detected 
any R° flux, even if only at their predicted 'relic' flux level, if the flux makes 
it to the detector. Then, the critical ingredient is the penetration depth l p of 
the R°; l p tends to be large because of a small energy transfer per collision. 
Using a geometric estimate of a(R°A) ~ ttA 2 ^ 3 (Fermi) 3 , Ref. j2(J estimates 
l p ~ 50 m x (m R o/l TeV). The result is that for m R o < 100 GeV, l p is too 
small for the flux to penetrate to any wimp detectors. Of course, l p would 
be further increased if the smaller a(R°A) estimate of were correct. On 
the other hand, the applicability of these limits is more marginal if the relic 
flux is as small as that obtained in Ref. p8| . Further, it is important to note 
that the effects of energy degradation are estimated to be small in Ref. |po| , 
whereas they are found to be sufficiently substantial in Ref. jl^| as to reduce 
rates below observable levels in all the cases considered there (assuming the 
fluxes are rescaled to the relic density level). 

(2) The capture and co-annihilation event rates depend upon the annihila- 
tion cross section and the fraction of the annihilation product energy appearing 
in the relevant channel, as well as whether or not the relevant density is halo 
or relic. For both types of observation, the conclusion seems to be that events 



would be likely to have been seen if the density is halo-like, but probably not 
if relic-like. 

(3) Since the R° can almost certainly be captured by oxygen ^| limits on 
exotic oxygen isotopes contained in water molecules in the ocean will place 
limits on the R°. The limits on such heavy water given in Ref. |2Cf| require 

(t>R° Si 0.1771ro( TeV) ^ 10 o M y \ cm~ 2 s _1 , where i acc is the accumulation time, 
with 100 Myr being an estimate for the age of the ocean. This flux limit can be 
compared with the fluxes computed earlier. One observes that if R 0, s are the 
dominant halo component, then we would require m R o > 100 TeV. However, 
we have already argued that the halo density should not be dominated by 
R 0, s. The R° flux estimated using the Ref. (2(J relic R° density prediction 
is well below the above limit. In Ref. (T^| the estimated relic R° density is 
even smaller, and the dominant source of R°'s impinging on the ocean is R° 
production in the atmosphere by cosmic rays. The crude estimate of the heavy 
isotope abundance arising from this source borders on the observed limit. 

In summary, it appears that we can escape the above limits for m R o < 
1 TeV so long as the halo density of i?°'s is similar to the (very small) relic 
density estimated in Ref. [Q. Alternatively, a relic density as large as that 
estimated in Ref. Q is allowed if energy degradation effects are much larger 
than allowed for in their wimp detection estimates or if m R o < 100 GeV. In 
either case, the i?°'s cannot be a significant source of dark matter. 



5 Beyond the MSSM 

Non-minimal models of low-energy supersymmetry can also be constructed. 
One approach is to add new structure beyond the Standard Model at the TeV 
scale or below. The supersymmetric extension of such a theory would be a 
non-minimal extension of the MSSM. Possible new structures include: (i) the 
supersymmetric generalization of the see-saw model of neutrino masses; (ii) 
an enlarged electroweak gauge group beyond SU(2)xU(l); (iii) the addition 
of new, possibly exotic, matter multiplets [e.g., a vector-like color triplet with 
electric charge |e; such states sometimes occur as low-energy remnants in Eg 
grand unification models]; (iv) the addition of low-energy SU(3)xSU(2)xU(l) 
singlets; and/or (v) the addition of one or more new families with massive 
neutrinos. 

A second approach is to retain the minimal particle content of the MSSM 
but remove the assumption of R-parity invariance pi|. The most general 



e We presume that the capture rate by protons is much smaller, if not zero. Otherwise, the 
much stronger limits on heavy hydrogen isotopes almost certainly rule out the R° LSP. 



explicit |^] R-parity-violating (RPV) theory involving the MSSM spectrum in- 
troduces many new parameters to both the supersymmetry-conserving and the 
supersymmetry-breaking sectors. Each new interaction term violates either B 
or L conservation. For example, new scalar-fermion Yukawa couplings would 
be: 

{^•L)pmnLpL m E n + (Xj^) pmn L p Q m D n + (Aij)p mTl U;D c m D c n , (14) 

where p, m, and n are generation indices, and gauge group indices are sup- 
pressed. In the notation above, Q, U c , D c , L, and E c respectively represent 
(u,cI)l, u\i d-L, {v, e~)i, and e c L and the corresponding superpartners. The 



Yukawa interactions are obtained from Eq. (14) by taking all possible combina- 
tions involving two fermions and one scalar superpartner. Note that the term 
in Eq. ( fli] ) proportional to Xb violates B, while the other two terms violate L. 

Additional R-parity violating effects arise from the mixing of L p and the 
Y = — 1 Higgs superfield Ha- This leads to new parameters \J p (p = 1,2,3) 
which are analogues of the supersymmetric Higgs mass parameter \i. Finally, 
there are soft-supersymmetry-breaking R-parity-violating A-terms and a B- 
term contributing to the scalar potential (i.e., the interaction of squarks, slep- 
tons and Higgs bosons): 

{0'L}pmnL p L m E n + (a,^)p mn LpQ m D n + {flB*)pmn 

U^D c m D c n + b' p L p H u , (15) 

where H u is the Y = 1 doublet of Higgs scalars. 

Phenomenological constraints on various low-energy L and B-violating pro- 
cesses yield limits on each of the coefficients (Az,) p „ m , (X' L ) pmn and (\B)pmn 
taken one at a time pl| . If more than one coefficient is simultaneously non-zero, 
then the limits are in general more complicated. All possible RPV terms can- 
not be simultaneously present and unsuppressed; otherwise the proton decay 
rate would be many orders of magnitude larger than the present experimental 
bound. One way to avoid proton decay is to impose B or L-invariance (either 
one alone would suffice). For example, models of "baryon-parity" possess a 
discrete symmetry that requires Xb = but allows non-zero L-violating terms 
in Eq. ph. ' 



6 The Higgs sector of low-energy supersymmetric models 

In the MSSM, the Higgs sector is a two-Higgs-doublet model with Higgs self- 
interactions constrained by supersymmetry p^ , |23] | . After absorbing unphysical 
phases into the definition of the Higgs fields one finds that the Higgs sector 

•''The phenomenology of models in which there is spontaneous R-parity violation from a 
sneutrino field vacuum expectation value is described elsewhere in these proceedings. 



is CP-conserving, so that tan/3 is a real parameter (conventionally chosen to 
be positive). The physical neutral Higgs scalars are CP-eigenstates. The five 
physical Higgs particles are those listed earlier. At tree level, tan/3 and one 
Higgs mass (usually chosen to be m A o) determine the tree- level Higgs-sector 
parameters. These include the other Higgs masses, an angle a [which measures 
the component of the original Y = ±1 Higgs doublet states in the physical 
CP-even neutral scalars], the Higgs boson self-couplings, and the Higgs boson 
couplings to particles of the Standard Model and their supcrpartners. 

When one-loop radiative corrections are incorporated, the Higgs masses 
and couplings depend on additional parameters of the supersymmetric model 
that enter via virtual loops. The impact of these corrections can be signif- 
icant. Most importantly, the tree-level MSSM-124 upper bound of m h o < 
mz\ cos2/3| < m z is increased to m h a < 125 — 130 GeV for m t — 175 GeV and 
a top-squark mass of m~ < 1 TeV. The charged Higgs mass is also constrained 
in the MSSM. At tree level, m 2 H ± = + m^ , which implies that charged 
Higgs bosons cannot be pair produced at LEP-2. Radiative corrections modify 
the tree-level prediction, but the corrections are typically small. 

In the parameter regime of m^o > 2mz, one finds that the couplings of 
the h° are nearly indistinguishable from those of the Standard Model Higgs 
boson. Moreover, the non-minimal Higgs states, H°, A , and are heavy 
and approximately degenerate in mass. This regime, called the decoupling 
limit, is rather generic and applies also to models with more general Higgs 
sectors. Present experimental bounds do not yet require the Higgs sector 
parameters to lie in the decoupling regime. In this regard, one of the most 
sensitive observables is the decay rate for b — > sj. Current data yields a rate 
that is within 2a of the Standard Model prediction. Applying these results to 
the MSSM, the observed rate (which is below the Standard Model prediction) 
would require a negative one-loop contribution of a light chargino and top- 
squark to cancel the positive loop contributions of the W ± and H ± . In the 
absence of light supersymmetric particle loops, one obtains a strong lower 
bound on the charged Higgs mass, since the latter contribution raises the 
predicted rate for b — > sj. In this case, one finds that m^o > 350 GeV, a 
result well into the decoupling regime. 

The MSSM Higgs mass bound quoted above does not in general apply to 
non-minimal supersymmetric extensions of the Standard Model. If additional 
Higgs singlet and/or triplet fields are introduced, then new Higgs self-coupling 
parameters appear, which are not significantly constrained by present data. 
These parameters can contribute to the light Higgs masses; the upper bound 
on these contributions depends on an extra assumption beyond the physics 
of the TeV-scale effective theory. For example, in the simplest non-minimal 



supersymmetric extension of the Standard Model (NMSSM), the addition of 
a Higgs singlet adds a new Higgs self-coupling parameter, A p4j. The mass 
of the lightest neutral Higgs boson can be raised arbitrarily by increasing the 
value of A. Under the assumption that all couplings stay perturbative up to the 
Planck scale, one finds in almost all cases that m^o < 150 GeV, independent 
of the details of the low-energy supersymmetric model [^5| . The NMSSM also 
permits a tree-level charged Higgs mass below mw- However, as in the MSSM, 
the charged Higgs boson becomes heavy and approximately degenerate in mass 
with A in the decoupling limit (where m^o 3> mz). 

Finally, in models of R-parity violating supersymmetry, the distinction be- 
tween the Higgs and matter multiplets is lost. Thus, R-parity violation permits 
the mixing of sleptons and Higgs bosons. The associated phenomenology will 
not be addressed here; see Ref. [M for further discussion. 



6. 1 Supersymmetric Higgs searches at future colliders 

We now consider the opportunities for detecting the Higgs bosons of R-parity 
conserving low-energy supersymmetry at future colliders. Over (nearly) the 
complete range MSSM Higgs sector parameter space, at least one of the MSSM 
Higgs bosons will be detectable at the LHC (if not at LEP-2 or the Tevatron), 
assuming that the machine runs at its design luminosity of 100 pb _1 per year, 
and under the assumption that the current detector design capabilities are 
achieved. An e + e~ collider (denoted eC) or a collider (denoted /j,C) 

would also guarantee detection of at least one Higgs boson for all of the MSSM 
Higgs sector parameter space (by extending the LEP-2 Higgs search) once the 
center-of-mass energy of the machine is above 300 GeV. 

However, complete coverage of the MSSM parameter space only means that 
at least one Higgs boson can be detected. For example, in the decoupling limit 
(where m^o > 2mz), h° will surely be discovered at LEP-2 or the LHC (and 
can be observed at the eC). But, detection of the heavier non-minimal Higgs 
states H°, A a , and/or H ± is not guaranteed. At the LHC, there is a region 
p7t , illustrated in Fig. |H in the (to^ , tan 0) parameter space characterized by 
tan 13 > 3 and m A o > 300 GeV | such that only the h° will be detectable. Q 
At the eC, H°A° and H + H~ pair production is not kinematically allowed 
if m^o > y/s/2. Moreover, the production rate for h°A° (although it may 

9 This is probably the only region of Fig. |s| that is consistent with the measured BF(b — * 37) 
given that r7ij"~ 1 TeV is assumed. 

h In the regime 1 <J tan /3 <J 3 and m A o 300 GeV, Ref. p7| asserts that A and H° can be 
discovered at the LHC via their tt decay mode. However, there has not been a full detector 
simulation and complete analysis of backgrounds to substantiate such a claim. 
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Figure 3: Discovery contours (5<r) in the parameter space of the minimal supersymmetric 
model for ATLAS+CMS at the LHC: L = 300 fb _1 . Two-loop/RGE-improved radiative 
corrections to the MSSM Higgs sector are included assuming m~ = 1 TeV and no squark 
mixing. 

be kinematically allowed) is suppressed in the decoupling limit. Thus, the 
non-minimal Higgs states are not detectable at the eC for m^o > \/~s/2. 

Recently, two important questions regarding the observation of the Higgs 
sector of a supersymmetric model have been addressed. 

• First, is discovery of one Higgs boson guaranteed if singlet Higgs fields 
are added to the MSSM? 

1. It has been demonstrated |^8| that at least one of the light CP-even 
Higgs bosons (^1,2,3) of a two-doublet plus one singlet NMSSM Higgs sector 
will be observed at the eC. This result follows from the fact that if the lightest 
(h±) has weak ZZ coupling then (one of) the heavier ones must be almost 
as light and have substantial ZZ coupling (and thus be discoverable in the 
e+e~ — > Z* — » Zh production mode). This result probably extends to models 



that contain several singlet Higgs fields. 

2. However, before the eC is built, only data from LEP-2, the Tevatron and 
the LHC will be available. One finds (2i| that discovery of at least one Higgs 
boson is no longer guaranteed if a singlet Higgs field is present in addition to 
the minimal two doublet fields. The no-discovery "holes" in parameter space 
are never terribly large, but are certainly not insignificant for moderate values 
of tan/3 ~ 3-10. 

In order to demonstrate point 2. above, one employs the same detection 
modes for the NMSSM as employed for Fig. | in the MSSM case: (1) Z* -> Zh 
at LEP-2; (2) Z* -> ha at LEP-2; (3) gg -> h -> 77 at the LHC; (4) gg -> 
ZZ* or ZZ A£ at LHC; (5) t -> H+b at the LHC; (6) gg -► 66/1, 66a -> 
b~b~T + T~ at the LHC; (7) gg h,a — > t + t~ at the LHC. Additional Higgs 
decay modes that could be considered, but not reliably, at the LHC include: 
(a) a — > Zh; (b) /i — ► oa; (c) /ij — > ^/ij; (d) a, /i — ► <t. In order to avoid 
having to consider these latter, we explicitly exclude any choice of NMSSM 
parameters for which (a)-(d) might be relevant. 

The free parameters specifying the Higgs sector of the NMSSM are: tan/3, 
the mass of the lightest CP-even Higgs mass eigenstate , the mass of the 
lightest CP-odd scalar m a (we consider only the parameter sub-region in which 
the 2nd CP-odd scalar is taken to be much heavier), the new superpotential 
coupling A appearing in W 3 XHdH u N , and three mixing angles ai.2.3 which 
parameterize the orthogonal matrix that diagonalizes the CP-even Higgs mass- 
squared matrix. One finds that A < 0.7 is required if A is to remain perturbative 
during evolution from mz to Mjj, which implies a tan /5-dependent upper limit 
on m,h 1 below ~ 140 GeV. 

As already noted, moderate tan (3 yields the worst case: for low tan j3 < 1.5, 
at least one of the CP-even Higgs is observable in the standard 77,4^,... 
modes; at high tan/3 > 10 — 11, the bbh and 66a modes (with h, a decaying to 
t + t~) are sufficiently enhanced to be observable. A typical unobservable re- 
gion of (ai, Q-ii 013) parameter space that arises at moderate tan/3 is illustrated 
in Fig. § 

3. At a collider, direct s-channel production — > hi is quite 
likely to be visible Only the (predicted) \i + [i~h\ coupling is needed; the 
h\ can be decoupled from ZZ, WW without affecting its detectability. If the 
hi has not been previously detected, then a scan search in the < 150 GeV 
region is required with AEbcam/^bcam 0.01%. The time required for this 
scan could run to many years, depending upon the luminosity of the muon- 
collider (/xC) and the size of the ji + \i~h\ coupling. 

• A second important question is: What is required in order to discover 
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Figure 4: For tan/3 = 5 and = 105 GeV, we display in three dimensional (01,0:2,0:3) 

parameter space the parameter regions searched (which lie within the surfaces shown), and 
the regions therein for which the remaining model parameters can be chosen so that no Higgs 
boson is observable (interior to the surfaces shown). 

the H°,A°, i? ± of the MSSM? 

Here, a region of particular concern is the decoupling regime of m^o > 
200 GeV, where finding evidence for a non-minimal Higgs sector will be the 
most problematical. This concern is most pronounced if tan [3 J> 3, i.e. such 
that the LHC cannot observe the H°,A° (see above). 

1. At the eC, the only means for detecting any of the heavier non- minimal 
Higgs states H°,A Q ,H ± when m^o > yfs/2 is to run the collider in the 77 
collider mode J3l],[32| . In this way, single 77 — > H° , ^4° production can be 
observed for m^o < 0.8y/s provided high integrated luminosity (e.g., L ~ 
200 fb _1 ) is accumulated. There is still no certainty that a 77 collider facility 
will be included in the final eC plans, nor any guarantee that the required 
luminosity can be achieved. 

2. At a collider, direct s-channel production ^+ — > H°, A may 



be observable p3| . Assuming no constraints on or knowledge of m^o , one would 
ideally wish to scan the entire region of m^o between y/s/2 (the pair production 
limit) and the maximum available yfs and be certain of seeing the H°, A signal 
(if present) for any tan j3 > 3 (the value below which the LHC will find the 
H , A ). For yfs ~ 500 GeV, this would require an integrated luminosity of 
L = 100 fb _1 (at modest beam energy resolution of 0.1%). However, current 
/LtC designs suggest that only L < 10 fb /yr will be achieved in this energy 
range. Thus, if we imagine devoting a period of several years to this scan, 
H°, A discovery at the /zC would only be certain if tan/3 > 5 — 6. As described 
shortly, the situation would be greatly improved if the h° has been detected 
and its properties measured at the eC (operating at yfs = 500 GeV) and/or 
the fiC (operating in the s-channel Higgs production mode). If only eC data 
is available, it would be possible to distinguish between the h° and hsM and 
at least roughly determine to^o so long as m^o is below 450 GeV. If we 
only have s-channel fj,C data the result is similar. If both types of data are 
available, a rough m^o determination would be possible up to m^o ~ 600 GeV. 
Precision measurements of the SM-like Higgs properties would thus allow a 
much narrower scan for H°,A° discovery; as a result, the H°,A° could be 
detected at the //,C with only one or two years running for essentially any 
tan > 1 - 2. 

3. At both the e + e~ and colliders, it is envisioned that the machine 

will be upgraded to increasingly larger yfs. Once yfs > 2m ^a, H°A° and 
H + H~ pair production will be observable. It has been shown that this will 
be true even if these Higgs bosons have substantial decays to supersymmetric 
particles 035). 

Returning to the MSSM Higgs sector, there are a number of interesting 
issues. 

• How precisely will we know m^o and how useful will it be? 

Certainly m^o will be very precisely measured: Am^o < 100 MeV or 
better is achievable at the eC and LHC; the fiC yields Am h o < 1 MeV via 
an s-channel scan. Precise knowledge of m h o has the potential to inaugu- 
rate a new precision game. At one-loop (neglecting mixing, etc.), Am^ = 

,24 mi 

g 2 —i In — |. For m/jO = 100 GeV, Am h o — ±100 MeV implies constraints of 

order Am t = ±170 MeV, Am~= ±2 GeV. (More generally and at two- loops, 
A and /i parameters enter.) This level of constraint on m~, . . . probably ex- 
ceeds that achievable experimentally by direct measurements. However, there 
is an important caveat: can theoretical predictions for m^o be made at the 
Am^o = ±100 MeV level? The currently available two- loop RGE-improved 
result is only reliable to a few GeV. 



• What are the prospects for and means for distinguishing the MSSM h° or 
other SM-like light Higgs bosons from the SM h$M via branching ratios 
and couplings, and what are the implications? 

We can only present a very brief summary. Further details can be found 
in fH and §7|. The Snowmass96 study |§ shows that LEP-2, TeV-33 and 
LHC data will not provide determinations of the Higgs couplings of sufficient 
accuracy to distinguish the h° from the hsM in the 'decoupling' regime of 
itia ^ 2mz- For example, for a SM-like Higgs mass ~ mz (the best case where 
data from all three machines will be available) the fundamental properties that 
can be determined (which is a limited subset) and their errors are: BF{bb) — 
±26%; (WWh SM ) 2 /(ZZ_hs M ) 2 — ±14%; (WWh SM ) 2 — ±20%; (ZZh SM ) 2 

- ±22%; {^h SM ) 2 /(bbhsM? - ±17%; SF( 77 ) - ±31%; (ggh SM ) 2 - 
±31%; (Uh SM ) 2 /(WWhsM) 2 — ±21%; (tth SM ) 2 — ±30%. At higher masses, 
LEP-2 data will not be available and TeV-33 statistics decline rapidly. Overall, 
without eC and/or fiG data, we will know that the Higgs is SM-like, but we 
will not be able to discriminate between different SM-like possibilities. 

As a next step, we imagine that we have data available from the eC as well 
as the LHC. This yields a dramatic improvement in our ability to determine 
the fundamental couplings of the h°. One finds |3(| that the best 'bet' for 
discriminating between SM-like Higgs bosons for mt < 130 GeV is via the 
ratio aBF(h — > cc) / aBF(h — > bb); for rrih > 130 GeV (as possibly relevant 
in the NMSSM), oBF(h -> WW*)/aBF(h -> bb) will be the most valuable 
ratio. The results from Snowmass96 use 'topological tagging' (e.g. primary, 
secondary, tertiary vertices for a b jet but only primary and secondary vertices 
for a c jet). For L = 200 fb _1 at -y/s = 500 GeV and m h < 130 GeV (and 
combining Zh, ZZ-iusion and VFW-fusion production processes), one obtains 
a statistical error in BF(cc) / BF(bb) of ~ ±7%. Theoretical uncertainty in 
m c (jn c ) and rrib(mb) (sum rules, lattice) and QCD running should reach the 
< 10% level in a few years, implying a net error of < 10%. The result is that one 
would be able to distinguish h° from Hsm at a > 2a level for m^o < 450 GeV 
for the typical m^o = mh SM = 110 GeV mass case. This is illustrated in Fig. |[ 
where it is seen that deviations are almost independent of tan f3 and are larger 
than - 20% for m A o < 450 GeV. 

For the NMSSM, the determination of the (ZZh) 2 coupling squared could 
be very valuable for checking the sum rule J2i(ZZhi) 2 — '^2 li (WWhi) 2 = 1, 
(where the (VVhi) 2 - V = W, Z - are defined relative to the SM- values). | A 
direct measurement (independent of the h branching ratios) is possible using 
e + e~ — > Zh with Z — > e + e~,/i + /i~ and e + e~ — > e + e~h (ZZ-fusion) produc- 



1 (ZZh) 2 is also crucial to getting many quantities in eC only analysis. 
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Figure 5: Constant value contours in (m ,tan 0) parameter space for the ratios 
[WW* /bb} h o /[WW* /bb]f lSM and [cc/bb\ h o /[cc/bb\f lSM . We assume "maximal mixing" in 
the squark sector and present results for the case of fixed m h o = mh SM = 110 GeV. The 
band extending out to large m A o at tan (3 ~ 2 is where m h o =110 GeV is theoretically dis- 
allowed in the case of maximal mixing. For no mixing, the vertical contours are essentially 
identical — only the size of the disallowed band changes. The same contours apply if bb is 
replaced by t+t~. 



tion and isolation of the ft as a peak in the missing mass recoiling against 
the £ + £~. Inclusion of the ZZ-fusion mode |58[ as well as the usual Zh as- 
sociated production mode, yields (ZZh) 2 errors in the rrih < 150 GeV range 
of order ±5% — ±7%. Determining (WWh) 2 is more involved. For exam- 
ple, for m/i < 140 GeV, one must measure a{yvh)BF(h — > bb), and divide 
by BF(h — * bb). One can determine BF(h — > bb), from the ratios of the 
type a(Zh)BF(h — > bb)/a(Zh) and the ZZ fusion analogue, with an accu- 
racy of about 5% for m/,. < 130 GeV (assuming a SM-like h). The result is 
- ±5% - ±6% error for (WWh) 2 when m h < 140 GeV. 

Determination of r^ ot (as needed needed to compute (Xh) 2 = BF(h — > 
V)r^ ot for X — bb, etc.) requires a lengthy indirect procedure for nih < 
130 GeV. One needs to determine BF(h — > 77) as well as possible and use 
77 collisions to determine T(h — > 77) and then compute r^ ot — gw^zp-^ ■ 
Unfortunately, BF(h — > 77) cannot be well measured at eC alone. The best 
results are obtained by using a complicated procedure involving LHC and eC 
data, as described in Ref. |p6| . The net result is an BF(h — > 77) error of 
< ±16% for rat 130 GeV. To measure T(hsM — > 77), the 77 collider 
is required. Current estimates, including systematics, suggest a FQism - > 
j~f)BF(h S M -> bb) error (for L = 50 fb^ 1 at the 77 collider) of 8%-10%. 
Given ±5% — ±6% error for BF(hsM — * bb) one finds a T(hsM — * 77) error 
of order ±12%. Putting all this together yields a r^ ot error of ~ ±18%; not 
very wonderful. A detailed summary of errors for branching ratios, couplings 
and the total width of a SM-like Higgs obtained by combining LHC, eC and 



77 collider data appears in Ref. |36 



Production of the h in the s-channel at a muon collider can provide even 
greater discrimination power between the h° and the Usm, provided 7^ 
m z fl33| . The crucial measurements are two: (1) The very tiny Higgs width: 
r^* = 1-10 MeV for a SM-like Higgs with m h < 140 GeV (i.e. mass 
as predicted for the h° of the MSSM). (2) (t(/x+ >~ -► h -> X) for X = 
T T , CC, bb,WW*,ZZ*. The accuracies estimated to be achievable during 
an optimized 3-point scan with total integrated luminosity of L = 0.4 ftV 1 (as 
currently projected for four years of operation with a beam energy resolution 
of R = 0.003%) in the case of m h o = 110 GeV are the following: t + t~ 
±8%; cc — ±19%; bb — ±3%; WW* — ±15%; ZZ* — ±190%; Tf* M - 
±16%. The important ratios for discriminating h from hsM, and their errors, 
are: ^ -» ±18%, ^= -» ±22%, W- -» ±15%, | -> ±20%. From 
Fig. ||, these errors imply that with L = 0.4 fb _1 at the fj,C one can distinguish 



■'Note that <x(/i+ /i — * h — * X) provides a determination of T(h — » )BF(h — » X) 

unless ct _ < rt ot . 



between the hsu and the h° at a > 2a level for m A a up to ~ 400 GcV, which 
is similar to the discrimination power of the eC operating at y/s = 500 GeV for 
L = 200 fb _1 . Note that since the r^ ot errors are big and since r^ ot is model- 
dependent, r^ ot is not so clearly useful as the above ratios. Still, deviations of 
r^o* horn are predicted to be substantial if m A a < 500 GeV. 

MSSM/SM l~(h— >/^u.) Ratio Contours 
m T0P = 1 75 GeV, m h = 110 GeV, No Mix 
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Figure 6: Constant value contours in (m A o , tan 0) parameter space for the ratio T(h° -> 
fi + H~)/T(hsM —* We assume "no mixing" in the squark sector and present results 

for the case of fixed m h o = m hsM = 110 GeV. For "maximal mixing", the vertical contours 
are essentially identical — only the size of the allowed parameter range is altered. Contours 
for F(h° -> bb)/T(h S M -» 66) are identical. 

Finally, let us consider what can be achieved if eC (or equivalent muon- 
collider) y/s — 500 GeV data (L = 200 fb _1 ), is combined with s-channel /xC 
data. The most powerful discriminator between the h° and Hsm turns out 
to be T(h — > This is because there are four statistically independent 

ratios that can be used to determine T(h — ► 

(1) [TjhsM -+ n + n-)BF(hs M -> bb)] M c , 

BF(/X SM - &&)eC 



1 ' BF(h SM -» WW*)eC 

1 j r(/ lSM - zz*) eC 

1 j r(/l SM - WW*) e C 

At = 110 GeV, the error for T(h — > /i + /i~) that results from combining 
all these determinations is ~ ±4%. (This error only increases to ~ ±5% if the 
fiC s-channel integrated luminosity is decreased from 0.4 fb^ 1 to 0.1 fb , for 
example.) Given that there is no systematic uncertainty in the prediction for 
r(h — > /i + /i~) (to m is precisely known), we see from Fig. ^ that T(h — > /i + /i~) 
provides a 3cr level of discrimination between the h$M and the h° all the way 
out to ra^o > 600 GeV. Similar results apply so long a m/, is not in the vicinity 
of mz ■ 

The importance of determining the SM-like Higgs properties with the 
greatest possible precision argues strongly for having both a eC and a /_tC, espe- 
cially since very substantial accumulated luminosity is needed at both machines 
in order to achieve good precision for a model-independent determination of 
all the Higgs couplings and its total width. 

• What can we learn from the branching ratios of the H°, A in the super- 
symmetric GUT context? 

Once the H°,A° and/or iJ ± have been detected (either in s-channel pro- 
duction at a fiC or in pair production at the eC), it will be possible to mea- 
sure their branching ratios to various final states with substantial precision. 
Refs. |^|J35| have demonstrated that such measurements, in combination with 
a determination of m^o and one of the gaugino masses, will allow determination 
of a large number of the MSSM parameters, thereby discriminating between 
different sets of initial Mjj boundary conditions. In particular, Ref. B4] shows 
that expected accuracies for supersymmetric mode decays in pair production 
will allow one to discriminate at a very high confidence level between different 
boundary condition models. Thus, Higgs pair production will be a powerful 
tool in determining the correct GUT-scale boundary conditions. 

Additional issues/questions that we do not have room to address in detail 
relate to a Higgs sector with CP-violation. If the Higgs sector is found to 
be CP-violating, then the MSSM is ruled out. However, the NMSSM and its 
extensions can have CP violation in the Higgs sector. Recent work in this area 
includes: 



• The development of a new sum rule for Higgs couplings and its use to 
extract limits from LEP-2 for a CP-violating general (non-MSSM) two 
Higgs doublet model (2HDM) or two-doublet plus one singlet (2D1S) 
Higgs sector |3{|. The sum rule can be used to show that LEP-2 data 
excludes the possibility that there can be two light Higgs bosons in the 
CP-violating 2HDM or three light Higgs bosons in the 2D1S model. 

• Development of optimal procedures for precision measurement of the 
couplings, including CP-violating couplings, for a neutral Higgs boson 
at a high energy e + e~ collider by combining iih and Zh production 
data |4C(| . One finds that truly excellent errors/limits on all couplings 
should be attainable. Prospects for determining the CP properties of 
a Higgs boson at the LHC using tth associated production are also not 
insubstantial |p~jj| . 

7 The lightest and next-to-lightest supersymmetric particles 

The phenomenology of low-energy supersymmetry depends crucially on the 
properties of the lightest supersymmetric particle (LSP). In R-parity conserv- 
ing low-energy supersymmetry, all Standard Model particles are R-even while 
their superpartners are R-odd. Thus, starting from an initial state involving 
ordinary (R-even) particles, it follows that supersymmetric particles must be 
produced in pairs. In general, these particles are highly unstable and decay 
quickly into lighter states. However, R-parity invariance also implies that the 
LSP is absolutely stable, and must eventually be produced at the end of a decay 
chain initiated by the decay of a heavy unstable supersymmetric particle. 

In order to be consistent with cosmological constraints, a stable LSP is 
almost certainly electrically and color neutral. Consequently (with the excep- 
tion of a gluino-LSP), the LSP in an R-parity-conserving theory is weakly- 
interacting in ordinary matter, i.e. it behaves like a stable heavy neutrino and 
will escape detectors without being directly observed. Thus, the canonical sig- 
nature for conventional R-parity-conserving supersymmetric theories is missing 
(transverse) energy, due to the escape of the LSP. 

As noted earlier, in mSUGRA models the LSP is likely to be the Xi ( as 
opposed to the gravitino) and it tends to be dominated by its U(l)-gaugino 
component (but, see below). A nearly pure U(l)-gaugino (sometimes called a 
bino) will be denoted by B. However, there are some regions of mSUGRA pa- 
rameter space where other possibilities for the LSP are realized. For example, 
there are regions of mSUGRA parameter space where the LSP is a chargino. 
These regions must be excluded since we reject the possibility of charged relic 
particles surviving the early universe. A sneutrino LSP is possible in mSUGRA 



models only if m~ < 80 GeV. In addition, the requirement that relic LSP's do 
not "overdose" the universe by contributing a mass density that is larger than 
the critical density of the universe rules out additional regions of mSUGRA 
parameter space. Requiring that the relic density of LSP's constitutes a signifi- 
cant part of the dark matter would impose further restrictions on the mSUGRA 
parameter space. 

In more general SUGRA models, the nature of the LSP need not be so 
constrained. One can envision a xJ-LSP which has an arbitrary mixture of 
gaugino and higgsino components^ by relaxing gaugino mass unification. A 
nearly pure wino LSP (denoted W) is the result if M2 < Mi. A nearly pure 
higgsino LSP (denoted by H) is possible in the region where the gaugino Ma- 
jorana masses satisfy Mi ~ M2 > The sneutrino can be a viable LSP (with 
no mass restriction as in mSUGRA models), although it is unlikely to be a 
major component of the dark matter. Finally, as described earlier, the gluino 
(more strictly speaking an R-hadron containing the gluino, e.g. R a = gg) can 
be the LSP §. 

In simple GMSB models, it was noted earlier that the mass of the gravitino 
lies in the eV-keV regime. Thus, in such models, the gravitino is the LSP, and 
the next-to-lightest supersymmetric particle (NLSP) also plays a crucial role 
in the phenomenology of supersymmetric particle production and decay. Note 
that unlike the LSP, the NLSP can be charged. In GMSB models, the most 
likely candidates for the NLSP are X? an d t^. ^ The NLSP will decay into 
its superpartner plus a gravitino [either Xi ~ * Ng 3 / 2 (N — 7, Z, or h°) or 
t r — > t ± 33/2], with a lifetime that is quite sensitive to the model parameters. 
In a small range of parameter space, it is possible that several of the next- 
to-lightest supersymmetric particles are sufficiently degenerate in mass such 
that each one behaves as the NLSP. In this case, these particles will be called 



co-NLSP's 42 1 



Different choices for the identity of the NLSP and its decay rate lead to a 
variety of distinctive supersymmetric phenomenologies 43 44|,|]|[l6 47, 48| , f4l| . 
These will be examined in detail later. 

A possible exception to the above gravitino LSP scenario is that discussed 
earlier. It is possible jl8| to construct a GMSB model in which the gluino 
is the LSP and the gravitino is much heavier (and, thus, phenomenologically 



fc In GMSB models, the tr tends to be the lightest slepton. It is lighter than cr and fiR 
because of the effect of the larger r-lepton-Higgs Yukawa coupling on the evolution of scalar 
masses from the messenger scale to the electroweak scale. Note that the same Yukawa 
coupling is responsible for tr—tl mixing, so that the light r mass eigenstate would not be a 
pure tr. Nevertheless, the tr component of the lighter r would dominate, so we ignore this 
distinction in the notation. 



irrelevant). The resulting phenomenology (to be discussed later) is not dissim- 
ilar to that associated with the SUGRA gluino-LSP case mentioned earlier in 
association with the O-II model studied in |J . 

Finally, we re-emphasize that restrictions on the SUGAR and GMSB 
model parameter spaces can arise from requiring that the relic LSP density 
not overdose the universe, but that it is premature to consider parameter re- 
strictions based upon requiring that the LSP be the primary component of 
dark matter; there are other possible sources of dark matter in all models. 



8 Classes of supersymmetric signals 

Due to a lack of knowledge of the origin and structure of the supersymmetry- 
breaking parameters, the predictions of low-energy supersymmetry depend on 
a plethora of unknown parameters. Many details of supersymmetric phe- 
nomenology are strongly dependent on the underlying assumptions of the 
model. Nevertheless, we can broadly classify supersymmetric signals at future 
colliders by considering the various theoretical approaches described earlier. 
In this section, we examine a variety of supersymmetric signatures, and in 
the next section we explore their consequences for experimentation at future 
colliders. 



8.1 Missing energy signatures 

In R-parity conserving low-energy supersymmetry, supersymmetric particles 
are produced in pairs. The subsequent decay of a heavy supersymmetric parti- 
cle generally proceeds via a multi-step decay chain , ending in the produc- 
tion of at least one supersymmetric particle which is weakly interacting and 
escapes the collider detector. Thus, supersymmetric particle production yields 
events that contain at least two escaping non-interacting particles, leading to 
a missing energy signature, f] (Only the missing transverse energy Ifo can be 
determined at hadron colliders.) 

In conventional SUGRA-based models, all supersymmetric events contain 
at least two LSP's which will escape the collider detector. This leads to the 
"smoking-gun" signature of low-energy supersymmetry: events with large miss- 
ing transverse energy in association with jets and/or leptons. However, there 

'We are aware of only two counterexamples. First, there are models in which the x? is 
the LSP and the lightest neutralino and chargino are nearly degenerate in mass. If the 
mass difference is <J 100 MeV, then x~t is long-lived and decays outside the detector fell pi. 
Second, there is the gluino-LSP scenario |],[l^j, the signatures for which will be elaborated 
upon later. 



are two unconventional approaches in which the smoking-gun signature is ab- 
sent. First, there is the model in which the Xi is the LSP but the lightest 
neutralino and chargino are nearly degenerate in mass. If the mass difference 
is < 100 MeV, then the Xi is long-lived and decays outside the detector |5j],[|. 
Second, there is the gluino-LSP scenario. Since a hadronic calorimeter mea- 
sures only the gluino kinetic energy, events in which the gluino (more precisely, 
the lightest R-hadron) is absorbed in or passes through the hadronic calorime- 
ter would also have substantial missing energy; see later discussion. However, 
there would be no jets arising from <? decays in such models. 

In conventional GMSB models, all supersymmetric events contain at least 
two NLSP's, and the resulting signature depends on the NLSP properties. Four 
physically distinct possible scenarios emerge: 

• (1) The NLSP is electrically and color neutral and long-lived, and decays 
outside of the detector to its associated Standard Model partner and the 
gravitino. 

• (2) The NLSP is the sneutrino and decays invisibly into ^53/2 either 
inside or outside the detector. 

In either of these two cases, the resulting missing-energy signal is then similar 
to that of the SUGRA-based models where Xi or v is the LSP. 

• (3) The NLSP is the Xi an d decays inside the detector to Ng 3 / 2 , where 
N = -f, Z or a neutral Higgs boson. 

In this case, the gravitino-LSP behaves like the neutralino or sneutrino LSP of 
the SUGRA models. But, the missing energy events in this case are character- 
ized by the associated production of (at least) two N's, one for each LSP. Q 
Note that if Xi is lighter than the Z (and ft ) then BR(x? -> 753/2) = 100%, 
and all supersymmetric production will result in missing energy events with at 
least two associated photons. 

• (4) The NLSP is a charged slepton (typically tr in GMSB models if 
m~ R < m~ ), which decays to the corresponding lepton partner and 
gravitino. 

If the decay is prompt, then one finds missing energy events with associated 
lcptons (taus). If the decay is not prompt, one observes a long-lived heavy 
semi-stable charged particle with no associated missing energy (prior to the 
decay of the NLSP). 

m If the decay of the NLSP is not prompt, it is possible to produce events in which one NLSP 
decays inside the detector and one NLSP decays outside of the detector. 



As noted earlier, there are GMSB scenarios in which there are several next- 
to-lightest supersymmetric particles (co-NLSP's) that are almost degenerate, 
in particular sufficiently degenerate that decays down to the true NLSP have a 
very long lifetime. The canonical example of this scenario is the case where the 
right-handed sleptons of the e and fi type have mass so close to that of the tr 
that they are effectively stable against chain decays down to the tr ■ Then any 
one of the three types of sleptons can emerge as the end product of a typical 
supersymmetric decay chain. The resulting signals can be linear combinations 
of some of the above scenarios. For additional details on the phenomenology 
of the co-NLSP's, see Ref. @. 

Let us return to the scenario in which a moderately massive gluino is 
the LSP. We have noted that this case can arise in both SUGRA models [§ 
and in non-conventional GMSB models |ll|. (It is important to note that in 
the GMSB case, the gravitino is expected to be much more massive than the 
gluino and will not enter into the phenomenology.) Supersymmetric events will 
terminate with the production of at least two gluinos that turn into the lightest 
R-hadron, most probably the R° = gg. Since a hadronic calorimeter measures 
only the R-hadron kinetic energy, events in which the R-hadron is absorbed 
in the hadronic calorimeter would also have substantial missing energy that 
would be an increasingly large fraction of the total energy as the gluino mass 
becomes large. However, since the g does not chain decay, the associated jet 
production in gg events would be very different from that of mSUGRA-based 
models. Experimental limits would need to be re-evaluated. 

In R-parity violating low-energy SUGRA models, the LSP is unstable. If 
the RPV-couplings are sufficiently weak, then the LSP will decay outside the 
detector, and the standard missing energy signal applies. If the LSP decays 
inside the detector, the phenomenology of RPV-models depends on the iden- 
tity of the LSP and the branching ratio of possible final state decay products. 
If the latter includes a neutrino, then the corresponding RPV supersymmetric 
events would result in missing energy (through neutrino emission) in associa- 
tion with hadron jets and/or leptons. However, other decay chains are possible 
depending on the relative strengths of Xl, X' l and Xb [see Eq. (Q)]. Other 
possibilities include decays into charged leptons in association with jets (with 
no neutrinos), and decays into purely hadronic final states. Clearly, these lat- 
ter events would contain little or no missing energy. If R-parity violation is 
present in GMSB models, the RPV decays of the NLSP can easily dominate 
over the NLSP decay to the gravitino. In this case, the phenomenology of the 
NLSP resembles that of the LSP of SUGRA-based RPV models. 



8.2 Lepton signatures 



Once supersymmetric particles are produced at colliders, they do not necessar- 
ily decay to the LSP (or NLSP) in one step. The resulting decay chains can be 
complex, with a number of steps from the initial decay to the final state p0| . 
Along the way, decays can produce real or virtual W's, Z's and sleptons, which 
then can produce leptons in their subsequent decays. Thus, many models yield 
large numbers of supersymmetric events characterized by one or more leptons 
in association with missing energy and with or without hadronic jets. 

One signature of particular note is events containing like-sign di-leptons 
p2| . The origin of such events is associated with the Majorana nature of the 
gaugino. For example, gg production, followed by gluino decay via 

9 mxt -> m^vxi (16) 



can result in like-sign leptons since the g decay leads with equal probability to 
either £ + or £~ . If the masses and mass differences are both substantial (which 
is typical in mSUGRA models, for example), like-sign di-lepton events will be 
characterized by fairly energetic jets and isolated leptons and by largely from 
the LSP's. Other like-sign di-lepton signatures can arise in a similar way from 
the decay chains initiated by the heavier neutralinos. 

Distinctive tri-lepton signals (H^j can result from xfx® ~ * (^ >± Xi)(^*Xi) 
when Z* — > £ + £~ . Such events have little hadronic activity (apart from initial 
state radiation of jets off the annihilating quarks at hadron colliders). These 
events can have a variety of interesting characteristics depending on the fate 
of the final state neutralinos. 

In GMSB models with a charged slepton NLSP, the decay £ — > £gz/2 (if 
prompt) yields at least two leptons for every supersymmetric event in associ- 
ated with missing energy. In particular, in models with a tr NLSP, supersym- 
metric events will characteristically contain at least two r's. 

In RPV models, decays of the LSP (in SUGRA models) or NLSP (in 
GMSB models) mediated by RPV-interactions proportional to Al and X' L will 
also yield supersymmetric events containing charged leptons. However, if the 
only significant RPV-interaction is the one proportional to \' L , then such events 
would not contain missing energy (in contrast to the GMSB signature described 
above). 



8.3 b- quark signatures 



The phenomenology of gluinos and squarks depends critically on their relative 
masses. If gluinos are heavier, they will decay dominantly into qq, ^] while the 
squark can decay into quark plus chargino or neutralino. If squarks are heavier, 
the squarks will decay dominantly into quark plus gluino, while the gluino will 
decay into the three-body modes qq\ (where \ can be either a neutralino 
or chargino, depending on the charge of the final state quarks). A number 
of special cases can arise when the possible mass splitting among squarks of 
different flavors is taken into account. For example, in mSUGRA models it is 
often the case that the third generation squarks are lighter than the squarks 
of the first two generations. If the gluino is lighter than the latter but heavier 
than the former, then the only open gluino two-body decay mode could be 
66. In such a case, all gg events will result in at least four 6-quarks in the 
final state (in associated with the usual missing energy signal, if appropriate). 
The 1, 24 and 200 Snowmass96 point scenarios of Table || yield precisely this 
situation. The g decay chains in these three scenarios are: 

-. ~90%T 7" 99% — OzT 79% — Ouu i / + - + - - -\ 

1: g -> b L b -> \2 bb ~> x\bb+(e^e ,vv,qq) 

S % ~0lT , ,T 

— ► Xi°b + bb 

_ . ~ 85% T T 70% ~o.T 99% . 0-O.T 2 §% ~0.T,T 

24 : g -> b L b -> x¥> h -> h u X ibb -> Xibbbb 



69% ~0~0~Oi.T 

XiXiXibb 



„„„ ~ 99% r T 100% ~ ,T 

200 : g — > b L 6 -> XiOO 



More generally, due to the flavor independence of the strong interactions, one 
expects three-body gluino decays into 6-quarks in ~ 20% of all gluino decays. [] 
Additional 6-quarks can arise from both top-quark and top-squark decays, and 
from neutral Higgs bosons produced somewhere in the chain decays p^ . 

These observations suggest that many supersymmetric events at hadron 
colliders will be characterized by 6-jets in association with missing energy ||,|35| . 



"In this section, we employ the notation qq to mean either qq or qq. 

"Although one top-squark mass-eigenstatc (ti) is typically lighter than b in models, the 
heavy top-quark mass may result in a kinematically forbidden gluino decay mode into tt\ . 
p Here we assume the approximate degeneracy of the first two generations of squarks, as 
suggested from the absence of FCNC decays. In many models, the 6-squarks tend to be of 
similar mass or lighter than the squarks of the first two generations. 



8.4 Signatures involving photons 

In mSUGRA models, most supersymmetric events do not contain isolated en- 
ergetic photons. However, some special regions of low-energy supersymmetric 
parameter space do exist in which final state photons are probable in the decay 
chains of supersymmetric particles. For example, if one relaxes the condition of 
gaugino mass unification [Eq. (0)], and chooses Mi ~ M2, then the branching 
ratio for X2 ~ > Xi7 can be significant [ p6[ . In one particular model of this type 
0, the x?-LSP is dominantly higgsino, while x.2 i s dominantly gaugino. Thus, 
many supersymmetric decay chains end in the production of X21 which then 
decays to Xxl- I n this picture, the pair production of supersymmetric particles 
often yields two photons plus associated missing energy. At LEP-2, one can 
also produce X1X2 which would then yield single photon events in association 
with large missing energy. 

In GMSB models with a x?-NLSP, all supersymmetric decay chains would 
end up with the production of Xi- in many models, the branching ratio for 
Xi — * 793/2 is significant (and could be as high as 100% if other possible two- 
body decay modes are not kinematically allowed). Assuming that Xi decays 
inside the collider detector, supersymmetric pair production would yield events 
with two photons in associated with large missing energy. The characteristics 
of these events differ in detail from those of the corresponding events expected 
in the model of Ref . || . 

8.5 Kinks and long-lived sparticles 

In most SUGRA-based models, all heavy sparticles decay promptly in the 
decay chain until the LSP is reached. The LSP is exactly stable and escapes 
the collider detector. However, exceptions are possible. In particular, if there 
is a supersymmetric particle that is just barely heavier than the LSP, then 
its (three-body) decay rate to the LSP will be significantly suppressed and it 
could be long lived. For example, in the models with Mi > M 2 |5l]|| implying 
m~± — m~ , the xt can be sufficiently long-lived to yield a detectable vertex, 
or perhaps even exit the detector. 

In GMSB models, the NLSP may be long-lived, depending on its mass 
and the scale of supersymmetry breaking, \/~F. The couplings of the NLSP to 
the helicity ±5 components of the gravitino are fixed by Eq. ([lj), as described 
earlier. For \/F ~ 100-10 4 TeV, this coupling is very weak, implying that 
all the supersymmetric particles other than the NLSP undergo chain decays 
down to the NLSP (the branching ratio for the direct decay to the gravitino is 
negligible). The NLSP is unstable and eventually decays to the gravitino. For 
example, in the case of the Xi-NLSP (which is dominated by its U(l)-gaugino 



component) one can use Eq. (Q) to obtain the decay rate T(xi — > 733/2) 



m-o cos 2 0w/167r_F 2 . It then follows that 



100 GeV \ ( VF 



For simplicity, assume that Xi ~> 733/2 is the dominant NLSP decay mode. 
If VF ~ 3000 TeV § then the decay length for the NLSP is cr - 100 m for 



m~ = 100 GeV; while Vf ~ 100 TeV implies a short but vertexable decay 

Xi 

length. 

A similar result is obtained in the case of a charged NLSP. Thus, if y/~F is 
sufficiently large, the charged NLSP will be semi-stable and may decay outside 
of the collider detector. 



8.6 Exotic super •symmetric signatures 

If R-parity is not conserved, supersymmetric phenomenology exhibits features 
that are quite distinct from those of the MSSM. Both AL = 1 and AL = 
2 phenomena are allowed (if L is violated), leading to neutrino masses and 
mixing, neutrinoless double beta decay, and sneutrino-antisneutrino mixing. 
Further, since the distinction between the Higgs and matter multiplets is lost, 
R-parity violation permits the mixing of sleptons and Higgs bosons, the mixing 
of neutrinos and neutralinos, and the mixing of charged leptons and charginos, 
leading to more complicated mass matrices and mass eigenstates than in the 
MSSM. 

Consequences for collider signatures are numerous. Most important, the 
LSP is no longer stable, which implies that not all supersymmetric decay chains 
must yield missing-energy events at colliders. All Xi decays contain visible 
particles: 

X?^ (3j), («('M, (&3,rtj) . (18) 

Thus, even e + e~ — » X1X1 P au " production becomes visible. A sneutrino de- 
caying via v — > vx\ also becomes visible. For Al 7^ 0, sneutrino resonance 
production in e + e~ J5^] or ti^ ' \i~ f5£| collisions becomes possible. For \' L 7^ 0, 
squarks can be regarded as leptoquarks since the following processes are al- 
lowed: e + u m — ► d n — > e + u m , vd m and e + d m — > u n — > e + d m . (Here, m and 
n are generation labels, so that c?2 = s, d$ — b, etc.) These processes have 



'Recall that larger values are probably forbidden in order to avoid overdosing the universe. 



received much attention during the past year as a possible explanation for the 
HERA high-Q 2 anomaly j5^]. The same term responsible for the processes 
displayed above could also generate purely hadronic decays for sleptons and 
sneutrinos: e.g., l~ — > u m d n and v p — > q m q n (q = u or d). If such decays 
were dominant, then the pair production of sleptons in e + e~ events would 



lead to hadronic four-jet events with jet pairs of equal mass 60 , a signature 



quite different from the missing energy signals expected in the MSSM. Alter- 
natively, At ^ could result in substantial branching fractions for t — » Iv 
and v — > £ + £~ decays. Sneutrino pair production would then yield events 
containing four charged leptons with two lepton pairs of equal mass. 



9 Supersymmetry searches at future colliders 

In this section, we consider the potential for discovering and studying low- 
energy supersymmetry at future colliders. Various supersymmetric signatures 
have been reviewed earlier, and we now apply these to supersymmetry searches 
at future hadron colliders. Ultimately, the goal of experimental studies of su- 
persymmetry is to measure as many of the MSSM-124 parameters (and any 
additional parameters that can arise in non-minimal extensions) as possible. In 
practice, a fully general analysis will be difficult, particularly during the initial 
supersymmetry discovery phase. Thus, we focus the discussion in this section 
on the expected phenomenology of supersymmetry at the various future facili- 
ties under a number of different model assumptions. Eventually, if candidates 
for supersymmetric phenomena are discovered, one would plan to utilize pre- 
cision experimental measurements to map out the supersymmetric parameter 
space and uncover the structure of the underlying supersymmetry-breaking. 



9.1 SUGRA-based models 

We begin with the phenomenology of mSUGRA. Of particular importance 
are the relative sizes of the different supersymmetric particle masses, which 
are predicted in terms of the mSUGRA parameters. An important generic 
feature of the resulting superpartner mass spectrum is that substantial phase 
space is available for most decays occurring in a given chain decay of a heavy 
supersymmetric particle. 

Extensive Monte Carlo studies have determined the region of mSUGRA 
parameter space for which direct discovery of supersymmetric particles at the 
Tevatron and the LHC will be possible (6^]. At the hadron colliders, the 
ultimate supersymmetric mass reach is most often determined by the searches 
for the strongly- interacting superpartners (squarks and gluinos), although at 
the Tevatron, the tri-lepton signature will be stronger in Run-II and at TeV-33 



if tan j3 is small. Cascade decays then lead to events with jets, missing energy, 
and/or various numbers of leptons. Gluino and squark masses up to about 400 
GeV (or the equivalent in terms of m~± , m~ ) can be probed at the upcoming 

Xi X2 

Tevatron Run-II; further improvements are projected at the proposed TeV- 
33 upgrade (6^], where supersymmetric masses up to about 600 GeV can be 
reached. The maximum reach at the LHC is attained by searching in the 
1£ + jets +i?T channel; one will be able to discover squarks and gluinos with 
masses up to several TeV |||. Some particularly important classes of events 
include: 

• pp — > gg — > jets +I$t and £ ± £ ± +jets-t^T, (the like-sign dilepton signal 
fl52[). For much of parameter space, the mass difference m ~~ m ^± can 



be determined from jet spectra end points [|52 55 64j], while Ti~± — ti~o 



can be determined from £ spectra end points in the like-sign channel 
64fl . An absolute scale for m~ can be estimated (within an accuracy 



55 



of roughly ±15%) by separating the like-sign events into two hemispheres 
corresponding to the two g's [ |52| , by a similar separation in the jets-b^T 
channel Ip3fl, or variations thereof B5l,64l. 



pp — » xf X2 ~^ (^ ± Xi)(-^*Xi)i which yields a trilepton +i?T final state 
when Z* — > £ + £~; m~ a — m~ is easily determined if enough events are 

X2 Xi 

available [ p3[ . 

-> « 2^ +^ T , detectable at the LHC for m~< 300 GeV @. 

Squarks will be pair produced and, for mo 3> TO1/2, would lead to 53 
events with two extra jets emerging from the primary q — > qg decays. 



The LHC provides significant opportunities for precision measurements of 
the mSUGRA parameters |55| . In general, one expects large samples of super- 
symmetric events with distinguishing features that allow an efficient separation 
from Standard Model backgrounds. The biggest challenge in analyzing these 
events may be in distinguishing one set of supersymmetric signals from an- 
other. Within the mSUGRA framework, the parameter space is small enough 
to permit the untangling of the various signals and allows one to extract the 
mSUGRA parameters with some precision. 

Important discovery modes at the eC include the following J65[ : 

• e+e- -> xtxi -> (qqXi or £vx°i) + (<?<7Xi or £z/£i); 

r The possible predominance of b jets has already been pointed out in Eq. (^). 



. e+e- -> (i+x? or 77x+) + (r*? or 

In both cases, the masses of the initially produced supersymmetric particles as 
well as the final state neutralinos and charginos will be well-measured. Here, 
one is able to make use of the energy spectra end points and beam energy 
constraints to make precision measurements of masses and determine the un- 
derlying supersymmetric parameters. Polarization of the beams is an essential 
tool that can be used to enhance signals while suppressing Standard Model 
backgrounds. Moreover, polarization can be employed to separate out various 
supersymmetric contributions in order to explore the inherent chiral structure 
of the interactions. 

The supersymmetric mass reach is limited by the center-of-mass energy 
of the eC. For example, if the scalar mass parameter mo is too large, squark 
and slepton pair production will be kinematically forbidden. To probe values 
of mo ~ 1-1.5 TeV requires a collider energy in the range of y/s > 2-3 TeV. 
It could be that such energies will be more easily achieved at a future 
collider. 

The strength of the lepton colliders lies in their ability to analyze super- 
symmetric signals and make precision measurements of observables. Ideally, 
one would like to measure the underlying supersymmetric parameters with- 
out prejudice. One could then test the mSUGRA assumptions, and study 
possible deviations. The most efficient way to carry out such a program is 
to first use low machine energies to study the light supersymmetric spectrum 
(lightest charginos and neutralinos and sleptons) , and make model- independent 
measurements of the associated underlying supersymmetric parameters. Once 
these parameters are ascertained, one can increase the machine energy and 
analyze with more confidence events with heavy supersymmetric particles de- 
caying via complex decay chains. Thus, the eC and LHC supersymmetric 
searches are complementary. 

Beyond mSUGRA, the MSSM parameter space becomes more complex. 
As discussed briefly earlier, it is possible to have non-universality among the 
scalar mass parameters (without generating phenomenologically unacceptable 
FCNC's) and/or non-universal gaugino masses. We discuss here only a few 
examples of non-universal gaugino mass scenarios. The first example is the 
Mi < M% possibility, for which the SU(2)-gaugino component is dominant in 
the lightest chargino and neutralino |5l]||. In this case, the Xi an d XT~ can be 
very degenerate, in which case the visible decay products in the xf ~ ► Xi • • • 
decays will be very soft and difficult to detect. ^ Consequences for chargino and 

s Of course, in the limit of extreme degeneracy, the X\ wn l be long lived and be detectable 
as a visible track in the vertex detector or electromagnetic calorimeter. 



neutralino detection in e + e~ and fi + /i~ collisions, including the importance of 
the e + e~ — > JX1X1 production channel, are discussed in Refs. plp[. There is 
also the possibility that m~ ~ m~± ~ m~ . The decay products in the g decay 
chain would then be very soft, and isolation of gg events would be much more 
difficult at hadron colliders than in the usual mSUGRA case. In particular, 
hard jets in association with missing energy would be much rarer, since they 
would only arise from initial state radiation. The corresponding reduction 
in supersymmetric parameter space coverage for Run-II at the Tevatron and 
TeV-33 is explored in Ref. §. 

As a second example, consider the case where the low-energy gaugino 
mass parameters satisfy M% ~ Mi . |^| If we also assume that tan f3 ~ 1 and 
\n\ < Mi, M2, then the lightest two neutralinos are nearly a pure photino 
and higgsino respectively, i.e., X2 — 7 an d Xi — H. For this choice of MSSM 
parameters, one finds that the rate for the one-loop decay ~* 7Xi dominates 
over all tree level decays of X2 an d BF(e — > ex°) ^ BF(e — > exi). Clearly, 
the resulting phenomenology differs substantially from mSUGRA expectations. 
This scenario was inspired by the CDF eejj event [ p4| . Suppose that the eejj 
event resulted from ee production, where e — ► ex\ — > cyxi- Then one would 
expect a number of other distinctive supersymmetric signals to be observable 
at LEP-2 (running at its maximal energy) and at Run-II of the Tevatron. In 
particular, at LEP-2 (s/s — 190 GeV) with L = 500 pb _1 of data one expects 
more than 50 U + X +I$t events and at least 3 77 + X +I$t events. At the 
Tevatron with L = 100 pb" 1 of data, one predicts more than 30, 2, 15, 4, 2 
and 2 events of the types U + X +tf T , 77 + X +tf T , £j + X +tf T , U~f + X +tf T , 
£jj + X+1$t, and ill + X +I$t , respectively. In the above signatures, X stands 
for additional leptons, photons, and/or jets. These signatures can also arise 
in GMSB models, although the kinematics of the various events can often be 
distinguished. 

As a final example, we return to the scenario in which the R = gg is 
the LSP. What would be the signatures at a hadronic accelerator? If m R o is 
not large, an initially energetic bare gluino will develop into a hadronic jet 
containing an R°. There will be little apparent missing energy since both 
the non-i?° components of the jet and the strongly-interacting R° itself would 
deposit energy in the hadronic calorimeter. Thus, at small m^o, current limits 
on g~g production at the Tevatron are no longer relevant due to the absence of 
missing energy. As m^jo increases, the R takes a larger and larger fraction of 

'We remind the reader that gaugino mass unification at the high-energy scale would predict 
M 2 ~ 2Mi [see Eq. (§)]. 

"To achieve such a small /i-parameter requires, e.g., some non-universality among scalar 
masses of the form m\, ^ mL, mL. 

d q I 



the gluino's initial energy and the energy of the R° as a fraction of m^jo that 
appears in the calorimeter becomes smaller and smaller. The maximum energy 
that can be deposited is the full kinetic energy, KE — m R a[l/y/\ — (3 2 — 1]. 
For full deposit of all the KE, the gluino would have to be effectively stopped 
in the detector. The rest mass energy would then effectively disappear, leading 
to substantial missing energy. However, from our earlier discussions, it must 
be remembered that the path length l p for the R° becomes large at large 
771 flu, so that an R° with mpo > 1 TeV is likely to simply pass through the 
detector leaving only a small amount of deposited energy, implying even larger 
missing energy. Whichever is the case, the Tevatron limits will again fail to 
apply, not because of an absence of missing energy, but rather because of 
the small amount of visible energy associated with the gluino-jet — the only 
visibly energetic jets would be those associated with initial state radiation. 
For intermediate (3 values, one could search for delayed energy deposits in the 
hadronic calorimeter (with an unusual profile associated with the fact that the 
R° would tend to deposit most of its KE at the stopping point). In either 
case, one could envision searching for the R using a high luminosity, high 
energy y/s J> several TeV beam dump experiment, in which mass spectrometer 
or chemical analysis of the dump could be used to search for the i?°-nuclear 
isotopes that would be formed. However, if juro is above 100 GeV or so, the 
predicted penetration or stopping length (as discussed earlier) would imply the 
need for a very large beam dump. 

Of course, if squarks are not too heavy, ^production will have substantial 
rate (if not at the Tevatron, then at the LHC) and each q will decay to qg. A 
heavy g-LSP will yield missing energy as described above and if m— — m~ is 
substantial there will also be significant associated jet activity. Tevatron limits 
employing the jets plusi^r signal would presumably apply. 

9.2 G MSB-based models 

The collider signals for GMSB models depend critically on the NLSP identity 
and its lifetime (or cquivalcntly, its decay length [which depends on the value of 
\[F as shown in Eq. (O)]). In this regard, it should be recalled that values of 
\[F as large as 10 3 — 10 4 TeV are certainly not disfavored in the GMSB context. 
Such large values of yf imply very long decay lengths for the NLSP. Thus, 
we examine the phenomenology of both promptly-decaying and longer-lived 
NLSP's. In the latter case, the number of decays where one or both NLSP's de- 
cay within a radial distance R is proportional to [1 — exp(— 2R/ct)] ~ 2R/(ct). 
For large cr, most decays would be non-prompt, with many occurring in the 
outer parts of the detector or completely outside the detector. To maximize 



sensitivity to GMSB models and fully cover the (y/~F, A) parameter space, we 
must develop signals sensitive to decays that are delayed, but not necessar- 
ily so delayed as to be beyond current detector coverage and/or specialized 
extensions of current detectors. 

In the discussion below, we focus on various cases, where the NLSP is 
a neutralino dominated by its U(l)-gaugino (bino) or Higgsino components, 
and where the NLSP is the lightest charged slepton (usually the tr) . We first 
address the case of prompt decays, and then indicate the appropriate strategies 
for the case of a longer-lived NLSP. 

• Promptly-decaying NLSP: — B 

We focus on the production of the neutralinos, charginos, and sleptons 
since these are the lightest of the supersymmetric particles in the GMSB 
models. The possible decays of the NLSP in this case are: B — > 7^3/2 or 
B — > Zg 3 /2- The latter is only relevant for the case of a heavier NLSP, and 
is, in any case, suppressed by a factor of tan 2 9w', it will be ignored in the 
following discussion. 

At hadronic colliders, the X1X1 production rate is small, but rates for 
Xtxi -> W^W^x°iXi -> WWwW77-fc#r, ^r^r -> i + ('X°iX°i -> ^~77 + 
$t, IiJl -> e + ^X°iXi -> e + e~JJ +rfr, etc. will all be substantial. Im- 
plications for GMSB phenomenology at the Tevatron have been studied in 
Refs. |f3| . [l8|j45| . [l^ . |66| . As noted earlier, it is possible to envision GMSB param- 
eters such that the ee77 -\-J$t CDF event corresponds to slepton pair produc- 
tion followed by Xi -> 753/2 However, in this region of GMSB 
parameter space, other supersymmetric signals should be prevalent, such as 
vl£l — * ^77+9r an d vlvl — * The X2X1 an d X1X1 rates would also 
be significant and lead to X77 +J$ T with X = £ ± , £+£' ~, £+£-£' ± . Limits on 
these event rates from current CDF and DO data already eliminate much of 
the parameter space that could lead to the CDF ee77 event [ |67f . 

At LEP-2/eC, the rate for the simplest signal, e + e~ — > Xi93/2 — * 7 
is expected to be very small [ ft7| . A more robust channel is e + e~ — > X1X1 ~~ * 
77 +I$t with a (flat) spectrum of photon energies in the range jy/s(l — (3) < 
E-y < jy/s(l + [3). Ref. |68| points out that there might be a small excess of 
such events in the current LEP-2 data, but the level of background requires 
careful evaluation . 

• Promptly decaying NLSP: x§ ~ H 

The possible decays of the NLSP in this case are: 
depending on the Higgs masses. If the corresponding two-body decays are not 



kinematically possible, then three-body decays (where the corresponding Higgs 
state is virtual) may become relevant. However, in realistic cases, one expects 
to contain small but non-negligible gaugino components, in which case 
X? — * 53/27 would dominate all three-body decays. In what follows, we assume 
that the two-body decay H — > g^nh? is kinematically allowed and dominant. 
The supersymmetric signals that would emerge at both Tevatron/LHC and 
LEP-2/eC would then be 46 + X +I$t final states, where X represents the 
decay products emerging from the cascade chain decays of the more massive 
superparticles. Of course, at LEP-2/eC direct production of higgsino pairs, 
e + e _ — > HH (via virtual s-channel Z-exchange) would be possible in general, 
leading to pure 46 +Ifcr final states. 

• Promptly decaying NLSP: £r 

The dominant slepton decay modes are: — > £^3/2 and £^ — » l^Xi* ~ * 

e ± (I%i T Y -> ^ ± (^ ± ^ T )'g 3 / 2 - The x? will first decay to ££ L and l£ R , followed 
by the above decays. 

At both the Tevatron/LHC and LEP-2/eC typical pair production events 
will end with £ R l R — > t + £~ +$t, generally in association with a variety of 
cascade chain decay products. The lepton energy spectrum will be flat in the 
£r£r center of mass. Of course, pure £r£r production is possible at LEP/eC 
and the £r£r center of mass would be the same as the e + e~ center of mass. 
Other simple signals at LEP/eC would include £l£l — * 6£ +^t- 

As noted earlier, if a slepton is the NLSP, it is most likely to be the tr. 
If this state is sufficiently lighter than the en and JIr, then en — > er^r and 
V-R ~~ * fJ-TRj decays (via the B component of the mediating virtual neutralino) 
might dominate over the direct br — * eg 3 / 2 and JIr —>■ /A93/2 decays and all 
final states would cascade to r's. The relative importance of these different 
possible decays has been examined in Rcf. |^0| . A study of this scenario at 
LEP-2 has been performed in Ref. fl49[ . 

• Longer-lived NLSP: i R 

If the £r decays mainly before reaching the electromagnetic calorimeter, 
then one should look for a charged lepton that suddenly appears a finite dis- 
tance from the interaction region, with non-zero impact parameter as measured 
by either the vertex detector or the electromagnetic calorimeter. Leading up to 
this decay would be a heavily ionizing track with (3 < 1 (as could be measured 
if a magnetic field is present) . 

If the £r reaches the electromagnetic and hadronic calorimeters, then it 
behaves much like a heavy muon, presumably interacting in the muon chambers 



or exiting the detector if it does not decay first. Limits on such objects should 
be pursued. Limits are currently quoted by DO and CDF for pseudo-stable 
charged strongly produced particles (stable quarks, squarks, . .. ), but the 
rates predicted in GMSB models for in pairs are very different. There will 
be many sources of £r production, including direct slepton pair production, 
and cascade decays resulting from the production of gluinos, squarks, and 
charginos. A detailed calculation is in progress Q. Based on the slepton 
pair cross section at the Tevatron, we estimate that a charged pseudo-stable 
£r can be ruled out with a mass up to about 80-100 GeV. Similar limits 
can probably be extracted from LEP-2 data. Including the other production 
mechanisms should yield a significant extension of the excluded mass range in 
many GMSB models. 

• Longer-lived NLSP: Xi 

This is a much more difficult case. As before, we assume that the dominant 
decay of the NLSP in this case is Xi ~> 753/2- Clearly, the sensitivity of 
detectors to delayed 7 appearance signals will be of great importance. As noted 
earlier, if \[F is so large that the Xi escapes the detector before decaying, then 
the corresponding missing energy signatures are the same as those occurring 
in SUGRA-based models. 

Consider first what is possible at the Tevatron without making use of 
the delayed 7 appearance signals. First, for any magnitude of v~F, the jets 
plus missing energy signal at the Tevatron will only be viable in the limited 
region of GMSB model parameter space |l(| characterized by the GMSB A 
parameter being < 30 TeV (< 40 TcV) in Run-II (at Tev-33); the standard 
tri-lepton signal will be viable for A < 65 TeV (75 TeV) |66| . The prompt-two- 
photon plus missing energy signal deriving from the many sources of prompt 
Xi — > 753/2 decays discussed earlier is also only visible in a limited portion 
of parameter space. The regions of viability for Run-II and at Tev-33 are 
illustrated in Fig. 0. Thus, even after combining the jets-plus-missing-energy, 
tri-lepton and prompt-two-photon signals, there will be a significant region of 
(V~F, A) parameter space that cannot be probed without using the delayed 
Xi -» 753/2 decays ©H). 

Thus, it is important that the collider detectors develop signals that are 
sensitive to the delayed Xi ~ * 753/2 decays. The ability to search for delayed- 
decay signals is rather critically dependent upon the detector design. The 
possible signals 0,^6) include the following: (a) searching for events where the 
delayed-decay photon is identified by a large (transverse) impact parameter as 
it passes into the electromagnetic calorimeter; (b) looking for isolated energy 
deposits (due to the 7 from a delayed Xi decay) in outer hadronic calorimeter 
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Figure 7: We present cross section contours in the [yF, A) parameter space for the (a) 
impact-parameter (ft), (b) outer-hadronic (OHC), (c) roof-array (RA), and (d) prompt- two- 
photon (27) signals. Contours are given at a = 0.16 fb (innermost), 2.5 fb (middle), and 
50 fb (outermost). Thus, the middle (outermost) contour corresponds to 5 events at Run- 
11/ L = 2 fb- 1 (Tev-33/L = 30 fb" 1 ). In this figure we have taken^? in = 50 GeV for all 
signals. 



cells, as possible with the DO detector; and (c) looking for delayed decays where 
the Xi decay occurs outside the main detector and the photon from the decay is 
observed in a scintillator (or similar device) array place at a substantial distance 
from the detector. The observed signal will always contain missing energy from 
the emitted 53/2 or un-decayed Xii which will be especially useful in reducing 
potential backgrounds. In all cases, timing information for the detected photon 
would be extremely valuable. Only if these (or closely related) techniques are 
employed can one detect GMSB supersymmetry in the bulk (and possibly 
preferred portion) of parameter space where neither \ F nor A is small. 

A first estimate of the utility of the above signals, is given in Fig. fj], taken 
from Ref. |6(| . In this figure we display cross section contours in the minimal 
GMSB model for the DO detector for the three types of signals described above. 
Cuts have been imposed that should eliminate backgrounds. (A DO collabo- 
ration study is needed for confirmation.) The cuts require a certain number 
of jets rij (defined with Ai? coa i = 0.5 and required to have pt > 25 GeV and 



\rj\ < 3.5) and one or more delayed photons (required to have E 1 > 15 GeV 
and be separated from jets by AR > 0.5) and a minimum amount of miss- 
ing transverse energy ^y- For all three delayed-decay signals (a), (b) and (c) 
presented in Fig. 0, we require rij > 3 and^r > 50 GeV. 

(a) The impact parameter (b) signal is defined by requiring that there 
be one or more Xi ~* 753/2 decays in which the impact parameter b of the 
photon (as measured using the pre-shower and electromagnetic calorimeter) is 
> 2 cm. 

(b) The outer hadronic calorimeter (OHC) signal is defined by requiring 
that there be one or more x? — ¥ 753/2 decays in central OHC cells. 

(c) The roof array (RA) signal is defined by requiring that there be one or 
more xi ~ * 753/2 decays in which the photon emerges at a vertical height above 
6.5 m (i.e. outside the detector, in particular, beyond the muon chambers) and 
passes through a 38 m x 28 m rectangular detection array centered vertically 
above the interaction point at a distance of 16.5 m (corresponding to the height 
of the roof of the DO detector hall). 

Assuming that the backgrounds to all these signals are negligible after our 
very strong cuts, 5 events will be adequate to claim discovery. Fig. |?] shows 
that in Run-II (the middle contours) the impact parameter signal will cover 
almost all of the ^f~F < 1500 TeV and A < 100 TeV portion of parameter space. 
The RA and OHC signals would provide important backup and confirmation 
for the impact parameter signal when A < 50 TeV. At Tev-33, the regions for 
which the b, OHC and RA signals are viable expand greatly. In combination, 
they provide sensitivity to GMSB supersymmetry throughout the preferred 
y/F < 3000 TeV, A < 150 TeV parameter region for the model examined. 
Clearly it is of great importance for the DO and CDF experiments to examine 
these signals more closely and develop reliable background estimates. 

Limits on delayed-decay photon signals from Run-I data (L ~ 100 pb~ ) 
require analysis of the OHC signal. A study is underway |7^] . The roof array 
was, of course, not present, and the impact parameter resolution was extremely 
poor due to the absence of a preshower. 

• Gluino-LSP 

The phenomenology for this variant jl8| of GMSB model would be much as 
described in the case of a SUGRA model with a gluino-LSP. The primary 
distinctions would arise from differences in the detailed mass spectra for the 
heavier sparticles. 



"Since the expected resolution in b is ~ 0.2 cm, this corresponds to a 10<r excursion from 
the primary interaction point. 



9.3 RPV models 



In R-parity violating models, the LSP is no longer stable. ^| The relevant 
signals depend upon the nature of the LSP decay. The phenomenology depends 
on which R-parity violating couplings [Eq. ([uj)] are present. For our discussion 
here, we assume that it is the Xi that is the LSP, although this is no longer 
required (by cosmology, etc.) if the LSP is not stable. 

• At the Tevatron and LHC 

If Xb 7^ 0, then Xi —> 3j (where j =jet). The large jet backgrounds imply 
that we would need to rely on the like-sign dilepton signal (which would still 
be viable despite the absence of missing energy in the events). In general, 
this signal turns out to be sufficient for supersymmetry discovery out to gluino 
masses somewhat above 1 TeV. However, if the leptons of the like-sign dilepton 
signal are very soft, then the discovery reach would be much reduced, [j] This 
is one of the few cases where one could miss discovering low-energy supersym- 
metry at the LHC. If Xl dominates Xi decays, Xi ~ > /i ± e T ^, e ± e T ^, and there 
would be many very distinctive multi-lepton signals. If X' L is dominant, then 
Xi — > tjj and again there would be distinctive multi-lepton signals. 

More generally, many normally invisible signals become visible. An impor- 
tant example is sneutrino pair production. Even if the dominant decay of the 
sneutrino is v — ► v\\ (which is likely if m~ > m~ ), a visible signal emerges 

from the xi decay as sketched above. Of course, for large enough or X' L 
the u's would have significant branching ratio for decay to lepton pairs or jet 
pairs, respectively. Indeed, such decays might dominate if m~ < "^~o- 

• At LEP-2, an e + e~ collider or a collider: 
The simplest process at the lepton colliders is: 

e+e- -> x°ix°i - mm, ggg, (tmm (19) 

Xb Xl X' l 

(or the collision analogue), where the relevant RPV coupling is indicated 

below the corresponding signal. Substantial rates for equally distinctive signals 
from production of more massive supersymmetric particles (including sneutrino 

™We assume that the gravitino is not relevant for phenomenology, as in SUGRA-based 
models. 

^Soft leptons would occur in models where m~± ~ m ~o (which requires non- universal 
gaugino masses Jsilpl). 



pair production) would also be present. All these processes (if kinematically 
allowed) should yield observable supersymmetric signals. Some limits from 
LEP data already exist Q . 

Of particular potential importance for non-zero is s-channel resonant 
production of a sneutrino in e + e _ |57| and | pq ] collisions. At a muon 

collider, this process is detectable down to quite small values of the appropriate 
Xl (due to the very excellent beam energy resolution possible at a muon col- 
lider) , and could be of great importance as a means of actually determining the 
R-parity- violating couplings. Indeed, for small R-parity-violating couplings, 
absolute measurements of the couplings through other processes are extremely 
difficult. This is because such a measurement would typically require the R- 
parity-violating effects to be competitive with an R-parity-conserving process 
of known interaction strength. (For example, R-parity-violating neutralino 
branching ratios constrain only ratios of the R-parity-violating couplings.) 
Since sneutrino pair production would have been observed at the LHC, an 
e + e~ collider and/or the (J. + fJ,~ collider, it would be easy to center on the 
sneutrino resonance in order to do the crucial sneutrino factory measurements. 

• At HERA: 

Squark production via R-parity violating couplings could be an explana- 
tion for the HERA anomaly at high x and large Q 2 (^9|. For example, if 
(A^,)ii3 ~ 0.04-0.1, a leptoquark signal from e + d^ — * £l ^ e + d would be 
detected if m~ < 220 GeV. 

• Delayed RPV decays: 

If the RPV coupling responsible for the decay of the LSP is sufficiently 
small, but not too small, the LSP decay might not be prompt but could still 
occur inside the detector or at least not far outside the detector. In this 
case, the same techniques discussed in the GMSB models for observing delayed 
Xi ~~ * 753/2 decays could prove very valuable. Delayed Xi decay OHC energy 
deposits might occur, leptons |^| from a delayed Xi decay might be observed (us- 
ing the electromagnetic calorimeter and its preshower) to have non-zero impact 
parameter, and the x? could decay outside the detector but before reaching 
the roof array. Especially at lower energy machines such as the Tevatron, a 
large extension in parameter space coverage would result from implementing 
these non-standard detection modes. Further, the relative abundance of events 

v If the RPV interaction leads to \i ~ > 3j, then the ability to measure the impact parameters 
of jets would be important, if possible with sufficient accuracy. 



obtained using these different techniques, especially if supplemented with tim- 
ing information, could allow a determination of the RPV coupling strength. In 
many scenarios, such determination would not be possible by any other means. 

Finally, we note that one should not rule out the possibility that RPV cou- 
plings could be present in GMSB models. Unless, the RPV couplings are rather 
small, the NLSP will decay via the RPV decay channels with a much higher 
rate than to channels containing the g 3 / 2 f^ll i and all the RPV phenomenology 
described earlier will apply. However, it is also possible that RPV and 53/2 
decays of the NLSP will take place at competitive rates. Delayed-decay signals 
could prove crucial to a full understanding of such a situation. 

10 Summary and Conclusions 

Low-energy supersymmetry remains the most elegant solution to the natural- 
ness and hierarchy problems, while providing a possible link to Planck scale 
physics and the unification of particle physics and gravity. Nevertheless, the 
origin of the soft supersymmetry breaking terms and the details of their struc- 
ture remain a mystery. There are many theoretical ideas, but we still cannot be 
certain which region of the MSSM-124 parameter space (or some non-minimal 
extension thereof) is the one favored by nature. The key theoretical break- 
throughs will surely require experimental guidance and input. 

Thus, we must rely on our experimental colleagues to uncover evidence 
for supersymmetry at present or future colliders. Canonical supersymmetric 
signatures at future colliders are well analyzed and understood. Much of the 
recent efforts, described in the previous section, have been directed at trying 
to develop strategies for precision measurements to prove the underlying su- 
persymmetric structure of the interactions and to distinguish among models. 
However, we are far from understanding all possible facets of MSSM-124 pa- 
rameter space (even restricted to those regions that are phenomenologically 
viable). Moreover, the phenomenology of non- minimal and alternative low- 
energy supersymmetric models (such as models with R-parity violation) and 
its consequences for collider physics have only recently begun to attract sig- 
nificant attention. The variety of possible non-minimal models of low-energy 
supersymmetry presents an additional challenge to experimenters who plan on 
searching for supersymmetry at future colliders. 

If supersymmetry is discovered, it will provide a gold mine of experimental 
signals and theoretical analyses. The many phenomenological manifestations 
and parameters of supersymmetry suggest that many years of experimental 
work will be required before it will be possible to determine the precise nature 
of supersymmetry-breaking and its implications for a more fundamental theory 



of particle interactions. 
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